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Abstract. We consider classical lattice models describing first-order phase transitions, and 
study the finite-size scaling of the magnetization and susceptibility. In order to model the 
effects of an actual surface in systems like small magnetic clusters, we consider models with 
free boundary conditions. 

For a field driven transition with two coexisting phases at the infinite volume transition 
point h = ht, we prove that the low temperature finite volume magnetization mfj-ee{L, h) per 
site in a cubic volume of size L'^ behaves like 

mfree(i, h) = + tanh(^^ L'^ {h - h^{L))) + 0(1/L), 

where h^{L) is the position of the maximum of the (finite volume) susceptibility and m± 
are the infinite volume magnetizations at h = ht + and h = ht — 0, respectively. We show 
that h^{L) is shifted by an amound proportional to 1/L with respect to the infinite volume 
transitions point ht provided the surface free energies of the two phases at the transition 
point are different. This should be compared with the shift for periodic boundary conditons, 
which for an asymmetric transition with two coexisting phases is proportional only to l/L^"^. 

One can consider also other definitions of finite volume transition points, as, for example, 
the position hu{L) of the maximum of the so called Binder cummulant U{j-ee{L,h). While 
it is again shifted by an amount proportional to 1/L with respect to the infinite volume 
transition point ht, its shift with respect to h^{L) is of the much smaller order We 
give explicit formulas for the proportionality factors, and show that, in the leading 
term, the relative shift is the same as that for periodic boundary conditions. 
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1. Introduction 



In the last twenty years, the study of finite size (FS) effects near first and second 
order phase transitions has gained increasing interest. While the study of FS effects for 
the second order phase transitions goes back to the work of Fisher and coworkers in the 
early seventies [FB72, FF69, Fi71], finite-size effects for first order phase transitions were 
first considered by Iniry [180] and, in the sequel, by Fisher and Berker [FB82], Blote 
and Nightingale [BN81], Binder and coworkers [Bi81, BL84, CLB86], Privman and Fisher 
[PF83], and others. 

Recently, these studies have been systematized in a rigorous framework by Borgs and 
Kotecky [BK90] (see also [BK92, BKM91]), and by Borgs and Imbrie [BI92a, BI92b, Bo92]. 
Their results cover both finite size effects in cubic volumes and long cylinders, both field 
and temperature driven transitions, but were always limited to periodic boundary condi- 
tions. While the periodic boundary conditions are natural for the description of computer 
experiments that are used to study the bulk properties of a system (note that periodic 
boundary conditions are used in these computer experiments because they minimize the 
unwanted finite size effects) they do not allow for the description of FS effects in actual 
physical systems like, e.g., small magnetic clusters, where surface effects are of major 
importance. 

In this paper we start a rigorous study of such surface effects. We consider spin systems 
in a finite box A = {l,-- - ,L}'^, imposing free or so called "weak" boundary conditions 
(see Section 2 below) instead of the periodic boundary conditions used in our previous 
work. 

In order to explain our main ideas, let us first review the FSS for a system in a periodic 
box [BK90, BK92, BKM91]. For a system describing the coexistence of two phases, say an 
Ising magnet at low temperatures, the partition function with periodic boundary conditions 
can be approximated by 



where Z± contain small perturbations of the ground state configurations cta = +1 and 
a A = —1, respectively. The error terms coming from the tunneling configurations can be 
bounded by 0{L'^e~^/^")e~^^^\ where f{h) is the free energy of the system and Lq is a 
constant of the order of the infinite volume correlation length. 

In the asymptotic (large volume) behavior of log Z± there should appear, in principle, 
volume, surface, and corner terms. A periodic box, however, has neither surface, 
nor edges or corners, and one obtains 



Zper{L, h) = Z. 



■+{L,h) + Z_{L,h), 



(1.1) 



Zper(L,/i) ^ e-M'^^'^" + e-^-W^"* = 2cosh 



f+{h) - f-{h) 




f+(h)+f_(h) 



2 



(1.2) 



where f+{h) and f-{h) are the (meta-stable) free energies of the phase plus and minus. 
Taylor expanding f±{h) around the transition point ht, and introducing the spontaneous 
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magnetizations m± of the phase plus and minus at ht, one obtains the FSS of the magne- 
tization mper{L, h) = L~'^dlogZpQj.{L, h)/dh in the form 

m++m_ , m+-m- f m+-m- ^^^^ 



mper(L, h) ^ ■■"^ ^ + ^ tanh (/, _ h,)L^ ^ . (1.3) 

It describes the rounding of the infinite volume transition in a region of width 

A/i ~ L-'^ (1.4) 

with a shift ht{L) — ht that vanishes in the approximation (1.3). A more accurate calcu- 
lation shows that, in fact, for a system describing the coexistence of two low temperature 
phases at the infinite volume transition point ht and with infinite volume susceptibilities 
X±i one has 

h^{L) -ht^ ^(X+-X-) r-2d ^ c)(^-3d) (15) 

(m_|_ — m_)'^ 

if h^{L) is defined as the position of the maximum of the susceptibility in the volume L'^. 

Turning to free boundary condition, we again expand logZ'-|-(L, h) into volume- surface- 
... -corner terms. This time, however, the volume A has a boundary, and the expansion 
yields 

-logZ±(L,/.) = f^\h)L'' + ft'\h)2dL''-' + OiL"-'), (1.6) 

where f±\h) — f±{h) are the (meta-stable) bulk free energies, while f^~^\h) are the 
(meta-stable) surface free energies of the phase plus and minus, respectively. As a conse- 
quence, (1.2) gets replaced by 



(1.7) 

At this point, one major difference with respect to (1.2) appears: while the free energies /+ 
and /_ are equal at the transition point ht , the surface free energies are typically different at 
ht (obviously, there are systems for which t+ := ^\ht) and r_ := ^\ht) are equal, 
as e.g. in the symmetric Ising model where r+ = r_ by symmetry, but for asymmetric 
first order transitions, this is typically not the case). The leading terms in the expansion 
around ht then lead to the formula 

,^ , , m_i_ -|- m_ m_i_ — m_ , / m_i_ — m_ , /^xx 
mt,UL,h) ^ + tanh i-^ {h - h^{L))L'^\ . (1.8) 



Here 
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which, for t_ 7^ t_(_, is now proportional to 1/L, while the width Ah of the transition is 
still proportional to L~'^. 

In fact, a formula of the form (1.8) has already been given in [PR90], with heuristic 
arguments very similar to those presented above. Here, our goal is twofold: first, we want 
to make the arguments leading to (1.8) rigorous, deriving at the same time precise error 
bounds on the subleading terms (in fact, our method allows to calculate in a systematic 
way the corrections to (1.8) in terms of an infinite asymptotic series in powers of l/L). 
Second, we want to generalize these results to a wider class of situations, including, in 
particular, the finite-size scaling of expectation values of arbitrary local observables. 

It will turn out that the more precise analysis of the subleading terms reveals an in- 
teresting fact: if one considers other standard definitions of the finite volume transition 
points, as e.g. the position hu{L) of the maximum of the so called Binder cummulant 
UireeiL, h), one finds that all of them are shifted, with respect to the infinite volume tran- 
sition point ht, by an amount proportional to 1/L. Their mutual shifts, however, are of 
the much smaller order l/L^'^, with proportionality factors that are the same as those for 
the corresponding shifts with periodic boundary conditions, see Section 2 for the precise 
statements. 

The finite-size scaling of local observables, on the other hand, will lead to the construc- 
tion of certain "meta-stable" states and their finite-volume analogues {■)±^, such 
that 

W£i , i±W+i^ ^ _ ,^(,)),.}. (1.10) 

Here A±{L) = (A)^'^* differ from the corresponding infinite volume expectation values 
A± = (A)'^ by an amount which is exponentially small in the distance dist(supp^, 9A), 
see Theorem 3.2 in Section 3.4 for the precise statement in the more general context of 
N phase coexistence. Note that the argument of the hyperbolic tangent in (1.10) is the 
same as in (1.9), and is independent of the particular choice for A. Thus the finite-size 
scaling of all local observables is synchronized in the sense that, after subtracting the 
"offset" ^ the functions {A)f^'^ asymptotically only differ by a constant factor, 

see Fig. 1. 

The organization of the paper is as follows: in the next section we present, in Theorem 
A, our main results for the finite size scaling of the magnetization and susceptibility in 
the context of a field driven transition with two coexisting phases. Section 3 is devoted to 
the contour representation of the models considered in Section 2, together with our main 
assumptions and results for a more abstract class of models describing the coexistence of 
N phases. We state two main theorems concerning the finite-size scaling: Theorem 3.1 
on partition functions and other ther mo dynamical quantities, and Theorem 3.2 on the 
finite-size scaling of local observables. In Section 4 we will construct suitable meta-stable 
free energies and prove Theorem 3.1, deferring the technical details to the appendices. In 
Section 5 we construct meta-stable states and prove the corresponding theorem. Theorem 
3.2. In Section 6 we prove the results stated in Section 2, using the abstract results 
formulated in Section 3. 
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Fig. 1. Finite size scaling of three different observables. 



2. Field driven transitions 



2.1. Definition of the model. 

In order to explain our main ideas, we consider an asymmetric version of the Ising 
model. Working on a finite lattice A = {1, . . . , L}'^, d > 2, we consider configurations 
cta:*!— >cr^e{ — 1,1} and the reduced Hamiltonian 



where J is the reduced coupling (containing a factor /3 = l/ksT), the first sum goes 
over nearest neighbor pairs (ij), while the third one is a finite range (i.e. ha — for 
diam^ < R, where R < oo) perturbation with translation invariant coupling constants 
ka e 1^- While the first two terms in (2.1) describe the standard Ising model, the third 
term is a perturbation that may break the +/— symmetry of the Ising model. We will 
assume that it is small in the sense that 




(2.1) 




where 6o > is a constant to be specified in Theorem A below. 
The partition function with free boundary conditions is 




(2.2) 
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The derivatives of its logarithm define the corresponding magnetization 

mtree{L, K) = L''^ ^ \og Ziree{L, h) (2.3) 

and the susceptibility 

XfreeiL, h) = —mhee{L, h). (2.4) 
The Binder cummulant, Uiree{L, h), is given as 

(r {M^)c _ 3((M - (M))^)^ - ((M - {M)r) 

where (•) denotes expectations with respect to the Gibbs measure corresponding to (2.1), 
(•)c denotes the corresponding truncated expectation values and M — ^j^y^ ^i- Note that 
Ui,ee{L, h) < 2/3 by the inequality (F^) > (F)^ (applied to F = (M - {M)f). 

2.2. Heuristic background, main ideas. 

For low temperatures (i.e. large J), the leading contributions to the partition function 
come from the constant ground state configurations cta = —1 and ua = +1- In this 
approximation, 

Zfree(i^, h) ^ e-^+(^''^) + e-^-^^'^\ (2.6) 

where 

F±(L,/i) = ^e±(i) (2.7) 



with the position dependent "ground state energies" 

1^ 
\A 



eot{i) = 2^ i^AjjY - ha, a = ±\. (2.8) 

ACA: 

In the same approximation, the magnetization mfree(-£') h) and susceptibility Xfree(-^) h) are 
given by 

mfree(L, h) - tanh( ~ ^+ ) (2.9) 

and 

Xfree(L,/.) - COSh"^ ( ^"^^ " ^+ ^-"^ ) . (2.10) 

Observing that ea(i) differs from the bulk value if i is in the vicinity of 9A, we 
expand E±{L, h) into a bulk term e±L^ plus boundary terms, 

F±(L, h) = e±ih)L'^ + e^l~^\h)2dL'^-'^ + 0(1'^-^). (2.11) 



Surface induced Finite Size Effects 



7 



While, still within the approximation by ground states, the bulk transition point ho is the 
value of h of which e+{h) = e-{h), the finite volume transition point ho{L) corresponds 
to the equality of E+{L, h) and E-{L, h). By (2.11), this leads to a shift 

hoiL)- ho = 0(1/1). (2.12) 

Notice that for periodic boundary conditions we get ho{L) = ho for zero temperature and, 
for nonvanishing temperatures, a shift ho{L) — ho proportional to l/L^'^ for periodic b.c. 
[BK90, BK92]. 

In order to make the above considerations rigorous, one has to take into account the 
excitations around the two ground states a\ = ±1. This is done in Section 3 and 4 and 
leads to a representation 

Zfree(i^, h) = (^e-^+(^'^) + e-^-^'^''^)) (1 + 0(L'^e-^/^o)), (2.13) 

where Lq is a constant of the order of the infinite volume correlation length and F±{L, h) 
have an asymptotic expansion similar to (2.11), namely 

F±{L,h) = f±{h)L'' ^ f'i-^\h)2dL''-^ ^0{L'-^), (2.14) 

where f±{h) are meta-stable free energies and fj^ ^\h) are (meta-stable) surface free 
energies. Once these results (see Theorem 3.1 in Section 3 for the precise statements) are 
proven, we obtain the desired finite-size scaling results by a rigorous version of the method 
presented in the introduction. 

2.3. Statements of results. 



In order to state our results in the form of a theorem, we introduce, for h ^ ht-, the free 
energy 

f{h) = f^^\h) = - lim L-'^logZfree(L,/i), 



the surface free energy 



/(rf-l)(/i) = - lim —r-TT 



logZi,,,{L,h)+L''f{h) 



., the corner free energy 



/W(/i) = - lim ^ logZfree(i:,/i) + i:"/(/i) + --- + 2'*-ML/W(/i) 



1 



L^cx) 2^^ 



d-1. 



(2.15a) 
(2.15b) 

(2.15c) 



as well as single phase magnetizations m± and surface free energies t± at the transitions 
point ht, 

m± = - —f{h) 



dh' 



ht±0 



r± = /('^-i)(/.,±0). 
We also recall that was defined as 



(2.16) 
(2.17) 



\K,\ 



A:OeA 



\A\ 
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Theorem A: Finite size scaling of m and x. 

Consider a perturbed Ising model with a perturbation of the form (2.1), with translation 

invariant coupling constants ka with range R < oo. Then there are constants Jq < oo and 
bo > such that, for \ \k\ \ < bo J and J > Jq , the following statements are true. Let 

AF{L) = f^'^-^\ht + 0)2dL^-^ + • • • + f^'^Hht + 0)2'^ 

- f^'^'^Hht - 0)2dL'^-^ - 0)2^ (2.18) 

and define h^{L) and hu{L) as the points where the susceptibility Xbee{,L, h) and the Binder 
cummulant U{ree{L,h) are maximal. Then^ 

mreeiL, h) = ^ + ^+ ~ tanh ^+ ~ - h^{L))L''^ + + ||^||)/L) 

(2.19) 

and 

X^reeiL, h) = (^^^^^) ^ COSh-^ ^ (^ - h^{L))L'^^ + 0((l + M)L'-') 

(2.20) 

provided \h - h^{L)\ < 0((1 + \\k\\)L-^). 

In addition, for AF{L) ^ 0, the shift h~^{L) obeys the bound 

= + r^^^^l Td^^ + ^^l/^))- (2.21a) 

777.-)- — TJl— J-J 

In the leading order, the shift of the point hu{L) with respect to ht is the same, 

hu{L) = ht+ ^^^i^(l + 0(1/L)). (2.21b) 

Remarks. 

i) If T+ ^ the equation (2.21a) (and similarly for (2.21b)) can be simplified to 

/, (L) = ht+ ^(1 + 0(1/L)) , 

yielding a shift ~ 1/L which is much larger then the width of the rounding, which, accord- 
ing to (2.19) and (2.20), is of the order l/L'^. 

ii) It it is interesting to consider the mutual shift h^{L) — hu{L). While both h^{L) — ht 
and hu{L) — ht are of the order 1/L, their mutual shift is actually much smaller, namely 

^(L) - h,(L) = 2 (^;:^:)3 + 0(i4t)- (2-22) 



■"^Here, and in the following, 0(L") stands for an error term which can be bounded by KL°', with a 
constant K that does not depend on h, J and k, as long as J > Jo and ||«;|| < bo J. 
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It is interesting to notice that, in the leading order l/L^*^, this mutual shift is exactly the 
same as the corresponding shift for periodic boundary conditions. 

iii) We stress that the condition \h — h^{L)\ < 0((1 + DL"-*^) is not a very serious 
restriction in our context, because the width of the transition in the volume L'^ is only 
proportional to L~^. In fact, in Section 6 we will close the gap left in Theorem A by 
showing that for \h — h^{L)\ > ^^^^ (1 + 1 l)-^""^, one has 

\miree{L,h)-m{h)\ < 0{1/L) (2.23) 

and 

\Xh,,{L,h)-x{h)\<0{l/L), (2.24) 

where m{h) and x(/i) are the infinite volume magnetization and susceptibility of the model 
(2.1). 

iv) Notice that, for periodic boundary conditions, it is possible to define finite size transi- 
tion points ht{L) with exponentially small shift, for example the point where mper(-£', h) = 
mper(2L, /i). Here, all these definitions lead to a shift ~ 1/L yielding no qualitative im- 
provement with respect to the point h^{L) or hu{L). 

v) In principle, the coefficients m±, t±, can be calculated up to arbitrary precision 
using standard series expansions, provided the microscopic Hamiltonian is known. On the 
other hand, the scaling (2.19), (2.20), and (2.21) would allow, in principle, to obtain the 
coefficients m_|_, m_ and the difference r_(_ — r_ from experimental measurements. 

vi) The general context considered in Section 3 allows to analyze the finite size scaling 
with more general boundary conditions then the free boundary conditions considered here, 
including, in particular, small applied boundary fields favoring one of the two phases near 
the boundary. In order to apply the techniques developed in this paper, it is necessary, 
however, to exclude boundary conditions which strongly favor one of the two phases. Such 
a condition is needed to ensure that the main contributions to the partition functions 
do in fact come from small perturbations of the two ground states a\ = ±1. For large 
boundary fields, the boundary may strongly favor one of the two phases. The leading 
contributions to the partition function then would include configurations which are in one 
phase near the boundary, and in the other one for the bulk. In such situations, wetting 
and roughening effects of the contour separating the boundary phase from the bulk phase 
would be important physical effects. We are not attempting to study these effects in the 
present paper. 
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3. General Setting and Main Theorem 

3.1. Contour Representation of the Ising Model. 

In this section we review the contour representation for the model (2.1). To make this 
subsection as simple as possible, and to have a concrete example at hand, we use for 
illustration the simplest symmetry breaking term, namely a perturbation of the form 

K (7i(7j(7k, (2.1') 

(iifc>cA 

where the sum goes over all triangles < ijk > made out of two nearest neighbor bonds < 
ij > and < jk >. See [PS75, 76] for the contour representation for the more general model 
(2.1). It will be convenient to introduce, in addition to the finite lattice A = {1, ■ ■ ■ , L}*^, 
the subset V — [|,L + ^]'^ ofM'^ which is obtained from A as the union of all closed unit 
cubes Cj with centers i E A. For a given configuration a\ C {—1, 1}^, we then introduce 
the set d as the boundary between the region V+ C V where = +1 and the region 
V- (ZV where ai = —1, and the contours Yi, - ■ ■ ,Yn corresponding to a a as the connected 
components of d. 

To be more precise, we define an elementary cube as a closed unit cube with a center in 
A (we sometimes use the symbol Cj to denote an elementary cube with center z e A), and 
introduce V± as the union of all closed elementary cubes Cj for which cTj = ±1, respectively. 
The set d is then defined as V+ fl V-, and the "ground state regions" V± are defined as 
V± \ d. With these definitions, the partition function with Hamiltonian (2.1') can be 
rewritten in the form 

d (TA 

where the second sum is over all configurations consistent with d. 

In order to specify the configuration a a, one has to decide which component of V \ d 
corresponds to CTj = +1 and which one to Ui = — 1. To this end, we introduce contours with 
labels. Given a configuration a a, the contours corresponding to a a are defined as pairs 
Y = (supp y, €({■)), where supp y is a connected component of d while a is an assignment of 
a label Q!(c) G {—1, +1} to each elementary cube that touches supp F^. It is chosen in such 
a way that a{ci) — ai. Note that the labels of contours corresponding to a configuration 
a A are matching in the sense that the labels a{c) are constants on every component of 
V\d. 

In fact, a set of contours {Yi, . . . , F„} corresponds to a configuration a a, if and only if 
i) supp Yi n supp Y^- = for i ^ j and 
ii) the labels of Yi, . . . , Y^ are matching. 



^In the language of [HKZ88], supp Y is called a (geometric) contour, while Y is called a labeled contour. 
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We call a set of contours obeying i) a set of non- overlapping contours and a set of contours 
obeying i) and ii) a set of non- overlapping contours with matching labels, or sometimes 
just a set of matching contours. 

In order to rewrite Ziy-ee{L, h) in terms of contours, we assign a weight p(Y) to each 
contour. This is done in such a way that 

n 

^-Hia,) ^ ^-E^iv^)^-E_iy-) -Q p^Yk). (3.1) 

fc=l 

Here Hia^) is the Hamiltonian (2.1'), Yi, . . . Y„ are the contours corresponding to ua and 

E±{V±)= J2 ^±W- (3-2) 

i€AnV± 

For the standard Ising model, p{Y) = e"'^'^', where |y| is the number of elementary (d—l)- 
dimensional faces in suppl". The third term in (2.1'), however, introduces corrections 
yielding a weight of the form p{Y) — e~'^^^^'^^^'^^^^\ As a consequence, 

|p(^^)| <e-^l^l withT = J-O(k). (3.3) 

Similar bounds hold for the derivatives \d'^p{Y)/dh'^\. 

With the help of (3.1), we rewrite the partition function Z{ree{L,h) as 

Z,UL,h)^ J2 e-^+(^+)e-^-(^-)np(n), (3.4) 

{Yi,...,Y„} k=l 

where the sum goes over all sets of matching contours in V . 

3.2. Assumptions for the General Model. 

In Section 3.3 below, we will state our main theorem. Theorem 3.1, from which we infer 
Theorem A of the preceding section. The setting of Theorem 3.1 is actually more general 
then what is needing for Theorem A and will include more general models. On one hand, 
we introduce contours in such a way that the notion of contours covers the Ising contours 
introduced above as well as thick Pirogov-Sinai contours [PS75, 76, Si82] constructed as 
unions of elementary cubes^. On the other hand, we also consider the situation of general 
N phase coexistence. 

As before, we consider the finite lattice A C Z*^, (i > 2, and the corresponding volume 
y C M''. We introduce the set C of elementary cells as the set of all elementary cubes in 
y, all closed d — 1 dimensional faces of these cubes, and all closed edges of these cubes. 



^The contours are introduced in such a way that the more general cases considered in [BW89, 90, 
HKZ88] are covered as well. 
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As usual, we define the boundary dW of a set C F as the set of all points x which have 
distance zero from both W and and as U dW. 

A contour in V is then a pair Y = (supp y, «(■)) where suppy is a connected union 
of elementary cells and «(■) is an assignment of a label a{c) from a finite set {1, . . . , N} 
to each elementary cube c in V \ supp Y which touches Y (by touching we mean that 
c n suppF 7^ while (c \ dc) fl suppF = 0). As before, we require that a is constant on 
each component C of F \ suppY, and say that a set {li, . . .Yn} of contours is a set of 
matching contours (or, more explicitly, a set of non-overlapping contours with matching 
labels) iff 

i) supp Yi n supp Yjy^^foriy^ j and 

ii) the labels of li, . . . are matching in the sense that they are constant on components 
of F \ (suppYi U • • • U suppYn)- 

In this way, each component C of V \ (supp Yi U ■ ■ ■ U supp Y^) has constant boundary 
conditions on dC\dV. The partition function of a statistical model with "weak" boundary 
conditions is then rewritten in terms of contours as 

n N 
{Yi,...,y„}fc=l m=l 

where the sum goes over sets of matching contours in V (including the empty set of 
contours), and Kn is the union of all components of F \ (supp U • • • U supp Y^) that have 
boundary condition m, and 

Em{Vm) = Yl ""rnic). (3.6) 

cCV^ 

We point out that the sum in (3.6) goes over all elementary cubes in the closure Vm of 
Vm, a convention which was chosen to ensure that all elementary cubes c with center in 
Vm are taken into account^. Note that by our definition of F as a closed subset of M^, the 
sum (3.5) contains contours that touch dV (in the sequel, we call these contours boundary 
contours, as well as contours that do not touch dV, ordinary contours). The contribution 
of the collection of empty contours to (3.5) is actually a sum of N terms, g-^m(^) 

In the equalities (3.5) and (3.6) we have introduced "contour weights" p{Y) e M and 
"ground state energies" Cm (c) e M that depend on a vector parameter h eU, where U is 
an open subset of M^. We assume that p{Y) and e^(c) are translation invariant as long as 
Y and c do not touch the boundary of V. More generally, we assume translation invariance 
along a, {d — k) dimensional face in dV as long as Y (or c) does not touch the {d — k) — 1 
dimensional boundary of this face. 

As usually, we have to assume the Peierls condition, together with several assumptions 
on the ground state energies em(c). Here, we assume that em{c) and p{Y) are functions 
of h obeying the following bounds: 

|p(y)| <e-^l^l-^°(^), (3.7) 

sum over elementary cubes c C Vm would exclude those elementary cubes c C Vm which touch one 
of the contours Yi,. . . ,Yn. 
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dh'' 



and 



< |A;|!(Co|r|)l'=le--l^l-^°(^), 



dh^ 



<a 



(3.8) 



(3.9) 



Here r > is a sufficiently large constant, \Y\ denotes the number of elementary cells in^ 

supp y , 

Eo{Y) = eo(c) with eo{c) = mmemic), (3.10) 

' m 

cCsupp Y 

A; is a multi-index k = {ka)a=i,...,i^ with 1 < |fc| < 6, \k\ = '^ka, and Cq is a constant 
independent of h and r. In addition, we assume that the difference between em{c) and the 
bulk term Cm is bounded, 

|em(c) - em\ < IT, (3.11) 

with a constant < 7 < 1 to be specified later. This condition is introduced to avoid a 
situation where free b.c. strongly favor certain phases n G {1, . . . , A^}. Note that 

|6m(c) Cm I ^ ll'^ll 

for the asymmetric Ising model (2.1). For this model, the condition (3.11) is therefore 
satisfied once < h^J for a suitable constant < 60 < 00. 

3.3. Main Theorem. 

In this section we state our main result for the general model introduced in the last 
section. It actually generalizes Theorem A presented in Sections 2 to a large class of 
models describing the coexistence of N phases. As in Section 2, the leading contribution 
to the partition function Z( V, h) is the sum 



m=l cCV ^ 



(3.12) 



Introducing \dkV\ as the joint A;-dimensional area of all A;-dimensional faces of V and e 
as solutions of equations 



(fc) 

m 



n=k ^ ^ 



(3.13) 



^Here, a fc-dimensional cell c in suppK is only counted if there is no (fc + 1) -dimensional cell c' in 
supp Y with c G c' . 
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whenever c is touching a fc-dimensional face of V and not touching its {k — 1) -dimensional 
boundary^, we rewrite 



J2 ^-(c) = e-l^l + ^m"'^ \dd-iV\ + • • • + e(^) \doV\. 



(3.14) 



To see that (3.13) imphes (3.14), just notice that a hypercube c touching a /c-dimensional 
face of V and not touching its {k — 1) -dimensional boundary is touching (^Z^) different 
n-dimensional faces of dV. Each of these faces is specified by choosing n — k directions 
among d — k directions orthogonal to the concerned A;-dimensional face. 
As usually we define the bulk free energy f{h) by 



f{h) = - \im^logZ{V,h). 

V^R'i \V\ 



(3.15) 



and the magnetization m{V,h) = {ma{V, h))ci=i,...,i/ by 

m{V,h) = ^J^log Z{V,h). 



(3.16) 



Theorem 3.1. There exist constants 6 > 0, 70 > and tq < 00 (where b and 70 de- 
pend on d and tq depends on d, N and the constant Cq introduced in (3.8) and (3.9) ), 
as well as meta-stable free energies fmih), surface free energies fm ^\h), ■ ■ ■ , edge free 
energies fm\h) and corner free energies fm\h), such that the following statements are 
true provided the effective decay constant f, 



T := r(l - 7/70) - To > 

(for the definition of t and 7 see (3.7), (3.8) and (3.11)). 
^) f{h) = min/^(/i). 

m 

ii) fm o-T^d fm , I = d — 1, . . . ,0, are 6 times differentiable functions of h. 



Hi) If \k\ < 6, then 



{fm (^m^ 



< e-^^ and 



d^ 



dh^ 



(0 



< e 



-hr 



dh^ 

where I — d — 1, . . . , 0. 
iv) Let 

FmiV, h) = fmih)\V\ + f^-'\h)\dd-iV\ + ■■■ + fj^\h)\doV\. 

Then 



(3.17) 



(3.18) 



dh'' 



N 



m—l 



< w 



\k\+l 



0(e 



-bfL\ 



maxe 

m 



-Fm(V,h) 



(3.19) 



^Note that due to translation invariance properties of em{c), the right hand side of this equation is 
constant for all such elementary cubes c. 
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provided < < 6. 
v) Let < < 5 and define Pg as 



N 



{V,h) 



■m=l 



-Fg{V,h) 



(3.20) 



Then 



d' 



dh'' 



m, 



N . 

Q=l ' ' ^ 



dFgjV.h) 
dhn 



< |V^|l'=IC>(e-''^^) 



(3.21) 



Here, as in the rest of this paper, 0{x) stands for a bound const a; where the constant 
depends only on d, N and the constant Cq introduced in (3.8) and (3.9). 

Theorem 3.1 is the main theorem of this paper. Its proof has three major parts: the 
geometric analysis of contours touching the boundary, a decomposition of Z{V, h) into pure 
phase partition functions, and the construction of meta-stable contour models allowing to 
prove the bounds (3.19) and (3.21). Deferring the technical details to the appendices, the 
main steps of this proof are presented in Section 4. 

3.4. FSS for Local Observables. 



In addition to the FSS of thermodynamic quantities like the magnetization or suscepti- 
bility, we want to study the FSS of local observables. In order to state our results in the 
context of the general models considered in Section 3.2, we introduce the following nota- 
tion. An observable A is a function which associates to each configuration contributing to 
(3.5) a real number A(Yi, • • • , Yn). Its expectation value in the volume V is defined as 



1 



Z{V,h) 



Z{A\V,h), 



where 



N 



z{A \v,h)= ^(^1' • • • ' n pi^k) n e-^-(^-) 

{yi,...,Y-„} 



(3.22) 



(3.23) 



fc=l m=l 



As in (3.5), the sum in (3.23) goes over sets of matching contours in y, and Vm is the 
union of all components of F \ (supp U • • • U supp Yn) that have the boundary condition 
m. 

An observable A is called a local observable, if there is a finite set of elementary cubes, 
denoted supp A in the sequel, such that A(Yi, • • • , Y^) does not depend on those contours 
Yi for which supp ^fl (supp Uint Yj) = 0, where Int Yi is the interior of Y"j (for the precise 
definition of Int Y"j see Section 4.1 below). 
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In most applications, local observables will be bounded, in the sense that the norm 

\\A\\= sup \A{Y^,--- ,Yr,)\ (3.24) 

{Y^,-,Y„} 

is finite. In addition, the observable will either not depend on the vector parameter h at 
all, or obey bounds of the form 

< |A;|!Cyi(Co|suppA|)l'=l, (3.25a) 

where Cq is the constant introduced in (3.8), Ca is a constant and /c is a multi-index of 
order < < 6. 

Here, we will allow for a slightly more general situation, requiring only that 

n 

< |A;|!C^(Co|suppyA|)'''l n^~^'^''~^°^^'^ ' (^-^^b) 

where supp Ya stands for the set supp A U supp Yi U • • • U supp Y^ , /c is a multi-index of 
the order < < 6, Co is the constant introduced in (3.8) and Ca is a constant that is 
finite^ for all h and r. Assuming this condition^ and the conditions introduced in Section 
3.2, we will be able to prove the following theorem. 

Theorem 3.2. There are "meta-stable expectation functionals" {•)v,q} Q = ' ' ' i ^ > such 
that the following statements are true provided the effective decay constant f := t(1 — 
7/70) ~ '!'o defined in Theorem 3.1. is positive and < \k\ < 6. 
i) For each local observable obeying the bounds (3.25a) or (3.25b), one has 

< C^e°(^)l'"PP^IO(e-^^-^), (3.26) 

where the probabilities Pq and the constant b are as in Theorem 3.1 and e = e""^/^. 
a) For each local observable obeying the bounds (3.25a) or (3.25b), the limits 

{A)l = lim^(A)^,, (3.27) 
exist as functions of h, and obey the bounds 

< C>(l)CA|suppA|l'=le°(^)l^"PP^I , (3.28) 



^While we assumed that the constant Cq is independent of h and r, we do not require that Ca is 
independent of h and r. 

*Note that (3.7), (3.8) and (3.25a) imply the bound (3.25b). 



A{Y^,.■■ ,Yn)llp{Yj 



^ 1 



dh^ 



N 



<7=1 



d" 



dh^ 
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where e = e 

in) For each local observable obeying the bounds (3.25a) or (3.25b), one has 



d 



< Ca I SUpp yl 1 1 e'^(^) I ""PP ^1 (e"''^ dist(supp A,dV) ^ 



where e = e "^Z^. 



Proof. The proof of Theorem 3.2 is given in Section 5. 
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4. Proof of Theorem 3.1 

The proof of Theorem 3.1 has three major parts: the geometric analysis of contours — 
in particular a bound of the form 

NgvilntY) < const |y I, 

where NQv{lB.t Y) denotes the number of elementary cubes in^ Int Y that touch the bound- 
ary dV of V, the decomposition of Z{V, h) into pure phase partition functions Zi{V, h), 
• • • , Zn{V, h), and the construction of suitable meta-stable free energies /i, • • • , fn- De- 
ferring the technical details to the appendices, we present the main steps in the following 
subsections. 

4.1. The geometry of contours. 

An important notion in the Pirogov-Sinai theory of contour models is the notion of the 
interior and exterior of a contour. For ordinary contours Y = (supp y, «(■)), one defines 
Int Y as the union of all finite components of \ supp Y and Int^ Y as the union of all 
components of Int Y which have the boundary condition m. Since ordinary contours do 
not touch the boundary dV of V, the set Ext Y = V \ (supp Y U Int Y) is a connected set 
and a{c) is constant for all cubes c in ExtF which touch suppF. We say that Y is an 
m-contour, if a{c) = m for these cubes. 

We now generalize these notions to boundary contours. To this end, we first introduce, 
for each corner k of the box F, an "octant" K{k). Namely, if k has components ki, . . .kd, 
with ki = 1/2 for i e I- and A;^ = L + 1/2 for z e /+, then 

K{k) := {x eR'^ \ Xi>l/2 for i e /_, < L + 1/2 for i e /+}. 

We then say: a contour Y is short iff there is a corner k such that suppY fl dV C dK{k). 
Otherwise Y is called long. Note that short contours may be ordinary contours or boundary 
contours, while long contours are always boundary contours. 

For a short contour Y , we then define Int Y as the union of all finite components of 
K{k)\ suppy, Intm Y as the union of all components of IntY which have the boundary 
condition m, Exty as F\ (suppY Ulnty) and V{Y) as suppY Ulnty. As before Exty is 
a connected set, and the notion of an m-contour is defined by the condition that a(c) — m 
for all cubes c in Exty that touch suppy. Note that these definitions are equivalent to 
the previous ones if the short contour Y is in fact an ordinary contour. Note also that the 
above definitions do not depend on the choice of the corner k if there are several corners 
k for which supp Y r\dV C K{k). 

For long contours, there is a priori no natural notion of an exterior or interior. We 
chose a convention that ensures that that the volume of a component Ci of Int Y cannot 
exceed the value L'^/2 if y is a long contour. Namely, if y is a long boundary contour. 



'We recall that we use the symbol W to denote the closure of a set W . 
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and Ci, . . . are the components of V \ supp Y, then the component Cj with the largest 
volume is called the exterior ExtY. If there are several such components C^^, . . .,Cii, we 
chose the first one in some arbitrary fixed order (for example the lexicographic order) as 
Ext Y. We then define Int y = F \ (supp Y U Ext Y) ,V{Y) = supp Y U Int Y, Int^ Y as the 
union of all components of Int Y which have the boundary condition m, and an m-contour 
y as a contour for which a{c) = m on all cubes c in Exty that touch suppy. 

The following three Lemmas state that the sets Ext Y and Int Y are defined in such a 
way that they have the main properties of an exterior and interior of the set suppy. They 
will be proven in Appendix B. 

The first of them expresses the fact that for two contours Yi and Y2, which do not touch 
each other, Yi together with it's interior is necessarily contained in one of the components 
of Exty2 Ulnty2. 

Lemma 4.1. Let Yi,Y2 be non-overlapping contours. Then the following statements are 
true. 

i) If supp Y2 C Extyi and suppyi C Extys, then V{Y2) C Extyi and V{Yi) C Exty2. 

ii) //suppyi C C2, where C2 is a component o/Inty2, then V{Yi) C C2. 
Hi) //suppyi C Inty2, then V{Yi) C lntY2. 

The next lemma expresses the fact that it is not possible that two contours which do 
not touch are both included in the interior of each other. 

Lemma 4.2. Let Yi and Y2 be non- overlapping contours. Then one and only one of the 
following three cases is true: 

i ) supp y2 C Ext yi and supp yi C Ext y2 , 

ii ) supp y2 C Ext yi and supp yi C Int y2 , 
Hi) suppy2 C Intyi and suppyi C Exty2. 

Definition 4.3. Let {yi,...,y„} be a set of non-overlapping contours. Then y^ e 
{Yi., . . . ,Yn} is called an internal contour iff there exists a contour Yi e {yi,...,y„} 
with supp Yk C Int Y^. Otherwise Yk is called an external contour. Finally, {yi, . . . , y„} is 
called a set of mutually external contours, if all contours in {yi, . . . ,Yn\ are external. 

The next Lemma will be used in Section 4.2 to conclude that all external contours of 
a given configuration contributing to (3.5) have the same external label. This observation 
will be an important ingredient in the decomposition of Z(y, h) into single phase partition 
functions Z^iV, h), and therefore in the proof of Theorem 3.1. 

Lemma 4.4. Let {Yi, . . . , Y^} be a set of non- overlapping contours in V, and let 



Then Ext is a connected component of V \ ljr=i ^^PP ^• 

Remark. Let Yq be a contour, and let Wq be one of the components of Int Yq. Then Lemma 
4.4 remains valid if V is replaced by Wq, as can be seen immediately from the proof in 
Appendix B. 



n 




(4.1) 
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While the preceding three Lemmas, even though tedious to prove, just express our 
intuitive notions about exteriors and interiors (in fact, our definitions were chosen in such 
a way that they do), the next Lemma is less obvious. In order to explain the need for it, we 
recall that the ground state energies ejn{c) may be different from the corresponding bulk 
term Cm- As a consequence, the boundary may favor an otherwise unstable phase. In the 
expansion about the leading contribution e~^"'^^^ to the single phase partition functions 
Zm{V, h), this will have the tendency to increase the weight of boundary contours which 
describe transitions into one of these "boundary favored" phases. In order to control the 
contributions coming from such contours (using the exponential decay e"^'^'), we need a 
bound of the form 

NgvilntY) < const 



where A^av(Int Y) denotes the number of elementary cubes in supp Y that touch the bound- 
ary dV of V. This is the main statement of the next Lemma. 

Lemma 4.5. Let Y be a contour in V, and let Wi, he the components o/IntF. 

Then 

NQv{\niY)<Ci\Y\, (4.2) 

n 

Y,\dWi\<C2\Y\, (4.3) 

i=l 

and 

\dV{Y)\<C^\Yl (4.4) 
where Ci = 2d{2^''^ + 1)7(2^/'^ - 1), Cs = Ci + 2d and C3 = Ca + 2d. 

The proof of this lemma relies on a lattice version of the isoperimetric inequality and 
is given in Appendix B. The proof of the required isoperimetric inequality is given in 
Appendix A. 

4.2. Decomposition of Z{V, h) into pure phase partition functions. 

The first step in the proof of Theorem 3.1 is the decomposition of Z{V, h) into N terms 
Zqiy, h), q = 1, . . . ,N , which are obtained as perturbations of the leading terms e~^«(^). 
We start with the observation that all external contours contributing to (3.5) touch the 
set Ext introduced in (4.1). Given that these contours are matching, we conclude that 
all external contours of a given configuration contributing to (3.5) have the same label. 
Therefore 

N 

Z{V,h) = J2z,iV,h), (4.5) 

q=l 

with 

n TV 

{yi,...,y„}/c=i m=i 



Surface induced Finite Size Effects 



21 



where the sum goes over sets of matching contours in V for which all external contours are 
Q'-contours. As before, Vm is the union of all components of F \ (suppYi U • • • U supp Y"n) 
that have boundary condition m, and Em{Vm) is defined in (3.6). 

More generally, let W he a, component of the interior Int Yq of some contour Yq in V, 
a set of the form (4.1), or a set obtained from a component Wq of an interior IntFo by a 
similar construction, 

n 

W = Wo\[j (Int Yi U supp Yi ) (4. 7a) 

i=l 

where {Yi, . . . , Y^} is a set of non-overlapping contours in Wq. We then define Zq{W, h) 
as 

n N 

{yi,...,y„}fe=i m=i 

where the sum goes over sets of matching contours in V for which all external contours 
are qf-contours with V{Y) = (suppY U Inty) C W. Here, Vm is now defined as the union 
of all components of W \ (supp li U • • • U supp Yn) that have boundary condition m. Note 
that the sum in (4.7b) contains no contours which surround the holes in W. Finally, given 
a volume W which is a disjoint union of volumes Wi, . . . , Wn of the form (4.7a), we define 
Zq{W, h) as the product of the partition functions Zq{Wi, h), i — 1, . . . ,n. 

Returning to (4.6), we derive a second expression for Zq{V,h), which eliminates the 
matching condition for the labels of Yi, . . . To this end we first sum over all sets 

{Yi, . . . , Yn} with a fixed collection of external contours. For each external contour Y this 
resummation produces a factor nm=i ^mi^^^m h). This yields the expression 



N 

{Yl ,. . .,K„ }ext fc=l 



PiYk) YI ^m(Intm>fe,/i) 



(4.8) 



where the sum runs over sets {Yi, . . . , l^jext of mutually external g-contours in V and 
Ext is the set defined in (4.1). Assuming that Zq{lntm Y^, h) ^ 0, we divide each Z^ by 
the corresponding Zg and multiply it back again in the form (4.7b). Iterating the same 
procedure on the terms Zq(Int^ Y^, /i), we eventually get 

n 

Z,(F,/.) = e-^«(^) YI n^'^(^'^)' (4-9) 
where the sum goes over set of non-overlapping contours which are all qf-contours, while 

i^,(r):=p(y).-.(.,|^^£^. (4,0) 

The equality (4.9) is the desired alternative expression for Zq{V,h) which contains no 
matching condition on contours. Assuming that the new contour activities Kq(Y) are 
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sufficiently small (for h in the transition region, this is actually the case, see Section 4.4), 
it also expresses the fact that e~^''^^^ is the leading contribution to Zq{V,h). 

Obviously, (4.9) can be generalized to volumes W of the form considered in (4.7). One 
obtains 

Z,(W^,/.) = e-^«W Yl fl^iiYk), (4.11) 

{yi,...,y-„}fe=i 

where the sum goes over sets of non-overlapping qf-contours Yi, . . . , y„ with V{Yi) C W. 
4.3. Truncated contour models. 

Given the decomposition (4.5) of Z{V,h) and the representation (4.9) for Zq{V,h), 
one might try to to obtain the FSS of Z(y, h) by a cluster expansion analysis of the 
partition functions ZqiV^ h). For such an analysis, one would need a bound of the form 
|-K'g(F)| < el^l with a sufficiently small constant e > 0. While it turns out that such a 
bound can be proven for stable phases it is false for unstable phases. 

In order to overcome this problem, we will construct truncated contour activities Kq{Y) 
and the corresponding partition functions 

n 

Z;(Ty,/i) = e-^^W WK^Y,) (4.12) 

{yi,...,y„}fc=i 

in such a way that: 

i) The truncated contour activities K'^iY) obey a bound 

|i^';(y)|<6i^i (4.13) 

for some small e > 0. 

ii) Z'qiW, h) = Zq{W, h) if the corresponding (infinite volume) free energy, fq = fq{h) is 
equal to / = min so that the truncated model is identical to the original model if 

fq^f (following [Si82] and [Z84], we caU these q "stable"). 

iii) The truncated contour activities, and the corresponding free energies, are smooth func- 
tions of the external fields h. 

Heuristically, the truncated model will be a model where contours corresponding to su- 
percritical droplets in the corresponding droplet model are suppressed with the help of a 
smoothed characteristic function. In the case of a two phase model, this idea could be 
implemented by defining 

K'^Y) = K+{Y)x{a\Y\ - (/+ - f-)\V{Y)\) 
K'_{Y) = K_{Y)x{a\Y\ - (/_ - f+)\V{Y)\), 

where x is a smoothed characteristic function and a is a constant of the order of r, for 
example a = t/2. While the presence of the characteristic function would not affect 
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the stable phases since f+ — f- > if + is stable (and f- — f+ > if — is stable), it 
would suppress contours immersed into an unstable phase + as soon as the volume term 
(/+ — is bigger then the decay term proportional to As a consequence, all 

contours contributing to the "meta-stable" partition function Z'^ obey a bound of the form 
(4.13) as desired. 

Unfortunately, the above definition of K'q{Y) is circular because it uses free energies 
fq that are defined as free energies of a model with activities K'^{Y). To overcome this 
problem, we will use the following inductive procedure. 

Assume that K'q{Y) has already been defined for all q and all contours Y with |y(y)| < 

n, n G N, and that it obeys a bound of the form (4.13). Introduce fq^ as the free energy 
of a contour model with activities 

[ otherwise. 

Consider then a contour Y with |y(y)| = n. Since |y(y)| < n for all contours Y in Inty, 
the truncated partition functions Z'^ilntmY, h) are well defined for all q and m. Their 
logarithm can be controlled by a convergent cluster expansion, and Z^(Intmi^, h) for 
all q and m. We therefore may define K'q(Y) for a qf-contours Y with |1^(F)| = n by 

if;(y) = 4(y)p(y)e-.(.,n^ii^, (4.i5a) 



with 



'.(Y) = n ^ ("1^1 - ifl^'^'^-f^-'^Win) ■ (4.15b) 



Here a is a constant that will be chosen later and x is a smoothed characteristic function. 
We assume that x has been defined in such a way that x is a function that obeys the 
conditions 

< x{x) < 1 , (4.16a) 
X{x) = if X <-l and x{x) = 1 if x > 1 , (4.16b) 

< -^x{x) < 1 , and (4.16c) 



< Co for aU k < 6, (4.16d) 



for some constant Cq. 

As the final element of the construction of K'^, we have to establish the bound (4.13) for 
contours Y with |F(y)| = n. We defer the proof, together with the proof of the following 
Lemma 4.6, to Appendix C. 
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We use fq — fq{h) to denote the free energy corresponding to the partition function 

Z'q{V,h), 

U = - hm -^logZ;(F,/i) , (4.17) 

and introduce / = f{h) and aq = aq{h) as 

f = minfm, (4.18a) 
aq = fq-f. (4.18b) 

Finally, we recall that for a volume W of the form (4.7a), \W\ denotes the euclidean volume 
of W, while for a contour Y or for the boundary dW of a volume W, \Y\ and \dW\ are used 
to denote the number of elementary cells, i.e. the number of elementary cubes, plaquettes, 
and bonds in Y and dW, respectively. 

Lemma 4.6. Assume that p(-) and eg(-) obey the conditions (3.7) and (3.11), and let 

e = e2+«e-(i-(i+2^i))^, (4.19a) 
2d 

a=-^{a-2). (4.19b) 

Then there exists a constant eo > (depending only on d and N) such that the following 
statements hold provided e < eo and a > 1 . 

i) The contour activities K'^{Y) are well defined for all Y and obey (4-13). 

ii) If aq\V{Y)\^''^ < a, then Xq{Y) = 1 and Kq(Y) = K'q(Y). 
Hi) If aq\W\^/'^ <a, then Zq{W,h) = Z'^{W,h). 

iv) For all volumes W of the form (4- 7a) one has 

\Zg{W,h)\ < e-/l^l)+o(e)|9w^l+7rWav.W^ (4_20) 

where NgviW) is the number of elementary cubes in W which touch dV. 

v) For W = V the bound (4-^0) can be sharpened to 

|^,(V^,/i)| <e-^l^le(i+^-)l^^lmax{e-^,e-(^^3)-V|av^||^ ^421) 

where C3 = Cz{d) is the constant defined in Lemma 4-5. 
Remarks. 

i) Due to the bound (4.13), the partition function Z'^iy, h) can be analyzed by a con- 
vergent cluster expansion. As a consequence, one can prove the usual volume, surface, 
corner asymptotics for its logarithm. Namely, using fq^\ fq^~^\ fq'^ to denote the 
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bulk, surface, comer free energies corresponding to Z'g{V,h), and introducing Fq{W) 
as 

F,{W) = J2 /.(c)' (4-22) 
cew 

where fq{c) = fqii c does not touch the boundary dV of and — in analogy to (3.13) 
— we have 

/^(^) = EClt)/-^' k = d-l,...,Q, (4.23) 

whenever c is touching a /c-dimensional face of V and not touching its {k — 1) -dimensional 
boundary, and as a result we get 

\\ogZ'q{V,h) + Fq{V) I < \V\0{{Ke)'^) (4.24) 

for some K < oo depending only on N and d. 

ii) It is interesting to present a heuristic derivation of the bound (4.21) in the approxi- 
mation of the droplet model. To this end, we recall that the diameter of a critical droplet 
is proportional to r/a^. Assume now that ciqL/r is small. Then the size of a critical 
droplet is larger then the system size, and Zq{V^ h) is a partition function describing small 
perturbations around a meta-stable ground state, with the weight 

Zq{V, h) ~ e-^-'l^l+o^l^^l) = e-^l^l+0(l^^l)e-««l^l. (4.25a) 

For large values of aqhjr^ on the other hand, supercritical droplets do fit into the volume 
V. As a consequence, the leading configuration contributing to Zq{V,h) contains a big 
contour (with an interior that is essential all of V) describing a transition from the unstable 
boundary condition q to a stable phase q with fq = f. We conclude that 

Zq{V, h) ~ e-^l^l+0(l^^l)e-OWI^^I (4.25b) 

if aqL/r is large. Except for the numerical value of the involved constants, the bound 
(4.21) exactly describes this behavior, 
iii) The fact that XqiX) suppresses supercritical droplets manifests itself in the fact that 

Xq{Y) = unless aq\V{Y)\ < (a + 1 + 0{e))\Y\ , (4.26) 

see Appendix C for the proof of (4.26). 

4.4. Bounds on Derivatives. 

We finally turn to the continuity properties of Zq and Z'q. As a finite sum of functions, 
Zq{V, h) is a function of h. The following lemma yields a bound on the derivatives of 

Zq{V,h). 
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Lemma 4.7. There is a constant K (depending on d, N and the constants introduced in 
(3.8), (3.9) and (4. 16)) such that the following statements are true provided e < eo and 
a>l. 

i) Zq{W,h) is a function of h and 



< |A;|!((Co + 0(e))|W^|)''^'e 



I'^l ^-f\V\^0{e)\dW\^'iTNgv{W) 



(4.27) 



for all multi-indices k of order 1 < |A;| < 6. 
a) K'g{Y) is a functions of h, and 



dh^ « 



< {Key 



(4.28) 



for all multi-indices k of order 1 < |A;| < 6. 
Hi) log Z'q{W,h) is a functions of h, and 



dh^ 



\og Z'{W,h) 



< + 0{e))\W\ 



(4.29) 



for all multi-indices k of order 1 < < 6. 
iv) For W = V (and 1 < \k\ < Q) , the hound (4-27) can be sharpened to 



d^ 
di? 



Z,{V,h) 



<\k\\ ((Co + Oie))\V\f^ e-^l^le(i+^-)l^^l 



{■ 



X max <! e "4,6 ^"^^^^ 



(4.30) 



Proof. The proof of this lemma is given in Appendix D. 
Remarks. 

i) For many models, including the perturbed Ising model introduced in Section 2, it 
is possible to prove a degeneracy removing condition. In the context of a model with A'" 
ground states and a driving parameter h e M^~^ (N = 2 for the perturbed Ising model), 
one considers the matrix 



^= {4ri^i-^N)] (4.31) 

q,i=l,...,N-l 



dhi 



and its inverse E ^ . One then proves that for some value Hq of h, all ground state energies 
are equal, and that obeys a bound 

||E-^||oo = max V |(E-^),J < const (4.32) 
q 
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in a neighborhood U oi ho, which does not depend on r. 
On the other hand, Sq = fq — Cq is a function of h with 



\U-eq\<0{e) 



and 



dh, 



< 0{e) . 



by Lemmas 4.6 and 4.7. As a consequence, the inverse of the matrix 



F 



_d_ 
dhj 



{fq - fN) 



(4.33) 



(4.34) 



(4.35) 



q,i=l,... ,N-1 



obeys a bound of the same form as E~^, with a slightly larger constant on the right hand 
side; combined with the inverse function theorem and the fact that /g(/io)~/iv(^o) ^ 0(e), 
one immediately obtains the existence of a point ht & U, \ht — ho\ < 0(e), for which all 
ttq are zero, i.e., all phases are stable. More generally, one may construct differentiable 
curves hq{t), starting at ht, on which only the phase q is unstable, surfaces hqq{t,s) 
on which phases q and q are unstable, etc. A possible parametrization of these curves, 
surfaces, etc., is given by am{hq{t)) = 5^q t, am{hqq{t, s)) = 5mq t + S^q s, 

ii) Due to Lemma 4.7 ii), the bound (4.24) can be generalized to the first six derivatives 
oflogZ'g{V,h). Namely, 



d 



{log Z'^{V,h)+Fq{V)) 



< \V\0{{Ke)^) 



(4.36) 



for all multi- indices k of order 1 < < 6. 
4.5. Proof of Theorem 3.1. 

In order to prove Theorem 3.1, we introduce the sets 

Q = {1,---,N} and S = {q e Q \ aqL <a}. 
Using the decomposition (4.5) together with Lemma 4.6 iii), we bound 



(4.37) 



d' 



N 



dh'' 



Z(F,/z)- J]e-^^(^) 

q=l 



qes 



+ 



q^S 



dh'' 



Zq{V,h) 



qiS 



dh'' 



(4.38) 



where k is an arbitrary multi- index of order < < 6. 
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Next, we observe that for 1 < \k\ < 6, 



F,{V) 



<0{l)\V\ 



(4.39) 



by the assumption (3.9) and the fact that FqiV) — EqiV) can be analyzed by a convergent 
expansion using Lemmas 4.6 and 4.7. 
For q & S, we then rewrite 



1 _ e^iiv)+iogZ'^{v,h) 



Using the bounds (4.24), (4.36) and (4.39), we obtain the foUowing bound on the first sum 
on the r.h.s. of (4.38), 



q€S 



Z'{V,h)-e 



<0((Ke)^)|F|l'=l+i^e-^^(^) < 



< 0((Ke)^)|y|l'=l+Vaxe-^<'(^^ < 0((Ke)^)|y|l'=l+^ maxe-^^^^^. (4.40) 



q£S 



In order to bound the last sum in (4.38), we observe that ior q ^ S one has 



dh'' 



< 0(l)|F|l^le-'^''(^) < C'(l)|F||fcle(^-+0(^))l^^le-/<'l^l 

= 0(l)|y|l'=le(^^+^('))l^^le-'*^l'^le--''l^l 

< 0(l)|V|l^le(27^+o(^»l^^le-«^l^l maxe-^«(^) 

~ q&Q 

< 0(l)|y|l'^le-(«/2'^-27-o(e))|aT/| niaxe-^''^^), 

qeQ 



(4.41) 



where we used the definition (4.37) of S and, in the last step, the fact that L '^\V\ 
{l/2d)\dV\. 

Finally, again for g ^ 5, we have 



d^ 



dh'' 



Zq{V,h) 



< 0(l)|V|l'=le(i+'^^)l^^lmax{e-«/4|^l,e-/4^3|ai-||e-/|v| < 

< 0(l)|F|l''le^^+^'^'^"'^'"^"/^'^''^/'^'^3})|5V'| jj^axe"'^''^^^ (4.42) 

qeQ 



by (4.21) and (4.30). 

Inserting the bounds (4.40) through (4.42) into (4.38), and observing that \dV\ > 2dL 
for all d > 2, we finally obtain the bound 



d' 



N 



dh^ 



q=l 



< 0(e-^/^°)|y|l'=l+Vaxe-^«(^^ (4.43) 

q&Q 
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where 

1 2d 
— := min{-log(Ke), a - 4ci7r - 0(e), a/4 -2d- M^r, —p^r - 2d- Ad'yr] 

2d 

= min{-log(Ke),a/4 -2d- M^r, -— r -2d- M^r}. (4.44) 

4C3 

Recalling the definitions (4.19) of a and e, together with the fact that C3 = Ci + M, we 
now rewrite 

- log(Ke) = r - a - (1 + 2Ci)7r - 0(1), (4.45a) 
2d 

a/A-2d-Ad^T= -— (a - (32d + 8Ci)7T) -O(l), (4.45b) 

4C3 



and 



_ 2d - 4d7T = ^ - (32rf + 8Ci)7t) - 0(1) 



= ^4^ (i " (16d + 4Ci)7T) - 0(1). (4.45c) 
Choosing a = ^ + {16d + SOi — 1/2)^t, we obtain 

- log(Ke) = ^ - (1/2 + 501 + 16d)7T - 0(1), (4.46a) 

a/4 -2d- M^T = f ^ - (1/2 + 16d + 5C'i)7t) - 0(1), (4.46b) 
4G3 V 2 / 

and 

= ^ " ^^/^ + + ^^^^^^) " ^^^^^ = 4^ - ^1 + + 10Oi)7r) - ro), 

(4.46c) 

where tq is a constant that depends on A*", d, and the constants introduced in (3.8), (3.9) 
and (4.16). 
Defining 

b = Hd) = j^. (4.47a) 

^° = ^"e** = i+32dVioc. - ^^-"""^ 

and 

f = t(1 - 7/70) - To, (4.47c) 
we obtain I/Lq = bf and hence the bound (3.19) of Theorem 3.1. 
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Observing that 

g = g2-f (i-(i+i6d+i0Ci)7) = £^£-5^^-"^/^°), (4.48) 



we note that the condition r > imphes the inequahty e < eo provided tq is chosen 
large enough. The condition a > 1, on the other hand, is trivial, since a — ^{a — 2) > 

It remains to prove statements iii) and v). While v) is a direct consequence of iv), iii) 

follows from the fact that (/m — em) and {fm — em) can be analyzed by a convergent 
cluster expansion involving the decay constant e. Observing that 0(e) < 0{e~'^) can be 
bounded by e"^"^, this proves iii). □ 
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5. Proof of Theorem 3.2 

5.1. Decomposition of Z{A\Vjh) into pure phase partition functions. 

The first step in the proof of Theorem 3.2 is the same as the first step in the proof of 
Theorem 3.1. Namely, we decompose Z{A\V,h) as 

N 

Z{A\V,h) = Y,ZM\V,h), (5.1) 

q=l 

with 

n TV 

Z,{A\V,h)^ J2 A{Y^,.-- ,Yr,)l[p{Yk)l[e-^-^''-\ (5.2) 

{Yi,...,Yn} k=l m=l 

Here the sum goes over sets of matching contours in V for which aU external contours are 
qf-contours. 

Next, we group all contours Yi for which V{Yi) fl supp^ ^ into a new contour Y^, 
and introduce the sets 

suppyA= U suppy, v{Ya)= U v{y), 

Y&Ya Y£Ya 

Int Ya = V{Ya) \ supp Ya , and Ext Ya = V\ V{Ya) , 

as well as 



supp Ya = supp Ya U supp A, V{Ya) = V{Ya) U supp A , 

Int(°) Ya = Int Ya \ supp A , and Ext^"^ Ya = Ext Ya \ supp A . 

As usual, Intm Ya is the union of all components of Int Ya which have boundary condition 
m, Int™ Ya = Int Ya n Vm, while IntJ^^^ Ya = Int^°) Ya n Vm- 

Recalling that A does only depend on those contours for which V{Y) fl supp A ^ 0, we 
then define 

n' AT 

p{Ya) = A{Yl, • • • , n P^^k) n e-'^'"(^-^(^^PP^\^"PP^^», (5.3) 

fe=l m=l 

where Y(, ■ ■ ■ , Y^, are the contours in Ya- Fixing now, for a moment, all contours Yi in 
(5.2) for which V{Yi) fl supp^ ^ 0, and resumming the rest, we obtain 

N 

Z,(A|y,/i) = 5^p(y^)Z,(Ext(°)y^,/i) H Zm{lnt^^^YA,h). (5.4) 

Ya m=l 
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Introducing 



K (Ya) = flfr.k'^^^^PP^-*) fr ^rnjlnt^^^YA.h) 



(5.6) 



we further rewrite (5.4) as 



Z,{A\V, h) = J2 i^(lA)e-^<'(™PP^-)Z,(Ext(°) Ya, /i)Z,(Int(°) Ya, h) . (5.6) 

Ya 



Z,(A|V,/^) = e-^^(^)5^i^,(y^) 5^ n^^(^'«)- 

{Yi,...,y-„}fc=l 



Using finally the representation (4.11) for Zq{Ex.t^^^ Ya, h) and Zg(Int*''^'' Ya, h), we get 

(5.7) 

Here the second sum goes over set of non-overlapping q contours Y\, . . . ,Yn, such that for 
all contours Yi, the set V{Yi) does not intersect the set suppY^. 

In order to make the connection to the standard Mayer expansion for polymer systems, 
we then introduce G{Ya, Yi, - ■ ■ , Y^) as the graph on the vertex set {0, 1, • • • , n\ which has 
an edge between two vertices i>l and j > 1, i 7^ j, whenever supp Yinsupp Yj 7^ 0, and an 
edge between the vertex and a vertex z 7^ whenever 1/(1^^) flsupp 7^ 0. Implementing 
the non-overlap constraint in (5.7) by a characteristic function (piYA, Yi, ■ ■ ■ , Y^) which is 
zero whenever the graph G has less then n-\-l components, the standard Mayer expansion 
for polymer systems (see, for example [Sei82]) then yields 



Z,{A\V,h) 
Z,{V,h) 



00 ^ 



Ya 



n=o {yi,...,y„} 



Y[K,{Yk) 



.{Ya,Y, 



1; • • • 



Yn) . (5.8) 



Here 0c (Y^, Yi, ■ ■ ■ , Y„) is a combinatoric factor defined in term of the connectivity prop- 
erties of the graph G{Ya, Yi, • • ■ , Y„), see [Sei82]. It vanishes if G{Ya, Yi, • ■ ■ , Y^) has more 
then one component. 



5.2. Truncated Expectation Values. 



In the context of Section 5.1, the expansion (5.8) is a formal power series in the activities 
K{Yi). In order to use this expansion, one has to prove its convergence. As in Section 4, 
it is useful to introduce truncated models. 

For a contour Y with V{Y) fl supp A = 0, we define K'^iY) as before, see (4.15a), while 
for Ya — {Yi, • • • , Y^}, where {Yi, ■ ■ ■ , Y^} is a set of contours with V{Y) fl supp A 7^ 
for all Y e Y^, we define 

k{ya) = p{YA)e-^^^^^ n n ' (^-q) 
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with Xg(^) cis in (4.15b). Given this definition, we introduce 



Ya {yi,...,y„}fc=i 



and 



Z'(V,h) 



which can again be expanded as 



oo ^ 

Ya n=0 ■ {yi,...,y„} 



n K^^k) 

.k=l 



(5.10) 



(5.11) 



(5.12) 



The following Lemma will allow us to prove absolute convergence of the expansion (5.12), 
which immediately yields the statements ii) through iv) of Theorem 3.2. 

Lemma 5.1. Let e, eo and a be as defined in Lemma 4-6, and assume that e < eo and 
a > 1 . Then the following statements are true: 

^) If 



thenKg{YA) = K'^{YA). 
ii) Let 



Then 



aq max \ViY)\^''^ <a . 



\Ya\= E |suppF| 

Y&Ya 



Hi) Let k he a multi-index of order 1 < < 6. Then 



< |suppyl|l'=ICAe^(^)l^"PP^I(Ke)l^^l , 



(5.13) 

(5.14) 
(5.15) 

(5.16) 



where K is a constant that depends only on d, N , and the constants introduced in (3.8), 
(3.9) and (4.16). 

Proof. The proof of Lemma 5.1 is given in Appendix E. 

Using standard estimates for polymer expansions, see for example [Sei82] , the bounds of 
Lemma 4.6 and Lemma 5.1 immediately imply the absolute convergence of the expansion 
(5.12), 



oo n 

\(^)vJ<J2\K(Ya)\J2-, E ll\Km\\MYA,Y^,---,Y^)\ 

< 0(l)CAe^('^l'"PP^I , 



(5.16a) 
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and similar bounds for the derivatives, in particular 



dh'' 



< 0(l)|suppA|I^IC^e^(^)l^"PP^I. 



(5.16b) 



Theorem 3.2 ii) through iv) then follows using standard arguments. 

5.3. Bounds on Zg{A\V,h). 

In conjunction with Lemma 4.6, Lemma 5.1 allows to analyze Zq{A\V, h)/Zg{V, h) pro- 
vided ttqL < a. In order to prove Theorem 3.2 in the case where ttqL > a for some of the 
phases q, we need an analogue of the bounds (3.21) and (4.30) for Zq{A\V, h). 

Lemma 5.2. Let e, cq and a be as defined in Lemma 4-6, le.t e = max{e, e""^"^/^}, and 
assume that e < eo and a > 1. Then the following statements are true: 



ii) Let k be a multi-index of order 1 < \k\ < 6. Then 

d^ 



dh' 



-ZqiA\V,h) 



< |A;|!((Co + 0(e))|F|)l'^ICAe°(^')l^"PP^I x 



Proof. The proof of Lemma 5.2 is given in Appendix E. 



5.4. Proof of Theorem 3.2. 

As pointed out before, the absolute convergence of the cluster expansion (5.12) immedi- 
ately implies the statements ii) through iv). In order to prove Theorem 3.2 i), we proceed 
as in the proof of Theorem 3.1, using the decomposition (5.1), Lemma 5.1 and Lemma 
5.2 instead of the decomposition (4.5), Theorem 4.6 and Theorem 4.7. Defining S as in 
Section 4.5, and observing that Zq{A\V, h) = Zq{A\V, h) ii q E S, we bound 



N 



dh'' 



q=l 

d'' 

q^S 



q€S 



q^S 



where k is an arbitrary multi- index of order < < 6. 



d^ 
dh^ 



(5.17) 
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Combining the bounds (4.40) and (5.16), and bounding terms of the form IsuppAjl'^l 
and by e'^^^^^, we get an estimate for the first sum on the right hand side of (5.17) 

by 



E 



Z'JA\V,h) - {A)^y^^e 



< CAe^^^^l^'^PP^I(Ke)^maxe-^<'(^) . (5.18) 



Here K is a constant that depends only on N, d and the constants introduced in (3.8), 
(3.9) and (4.16). The terms ior q ^ S are bound in a similar way, leading to 



d^ 
dhF 



< C^eO(^)l^"PP^Ie-(«/2^-2^-0(i))l^^l maxe-^''^^) (5.19) 

q&Q 



and 



d^ 



dh^ 



Z,{A\V,h) 



qeQ 



Inserting the bounds (5.18) through (5.20) into (5.17), and choosing a as in Section 4.5, 
we get e = e and 



d'' 

dhF 



N 



q=l 



< C^e^(')l'"PP^IO(e 



-^^•^)maxe-^<'^^) 
qeQ 



(5.21) 



where b and f are the constants introduced in (4.47). Together with Theorem 3.1 and 
the observation that a prefactor can be absorbed into the exponential decay term 

e-''^^, the bound (5.21) implies Theorem 3.2 i). □ 
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6. Proof of Theorem A 

Even though the statements of this Section have a generahzation (sometimes a very 
straightforward one) to the case of several phases, we will restrict ourselves to the situation 
where only 2 phases, plus and minus, come to play and the driving parameter is an external 
field h. However, we do not restrict ourselves to the model (2.1) (for which Theorem A is 
stated), but consider the two phase case in the general setting of Section 3. In particular, 
we have two ground state energies e± satisfying, for h in an interval U (containing the 
point ho for which e+{h) — e_(/i)), the (nondegeneracy) bounds 

0<a<-^{e_{h)-e+{h))<A (6.1) 

(Jjf b 

that imply the bounds 

Q<a<^{f_{h)-U{h))<A (6.2) 

on the free energies /±(^) from Theorem 3.1 (cf. also (4.17)). Actually, A = 2Cq according 
to the assumption (3.9). In the situation of Theorem A we have -^{e-{h) — e+(h)) = 2. 
Considering now the free energies^° 

d 



F±{L,h) = J2f±ih)\dkV\, (6.3) 

cf. (3.18), and their derivative 

M±{L,h) = J2^±\h)\dkVl (6.4) 



k=0 



d 



k=0 

where m^^\h) = —df^\h)/dh, and introducing 

AF{L,h) = F+{L,h) - F^{L,h) , (6.5a) 
AM{L,h) ^ M+{L,h) - M_{L,h) (6.5b) 

and 

M„(i.ft) = M±(M)+^Mi^, ^^^^^ 
the bounds (3.21) of Theorem 3.1 can be, for the two phase case, reformulated as 

[ ,^ /I 1 AM(L,/i) , , AF(L,/i) 



m{L,h) - [j^Mo{L,h) + ^ tanh( ^)j 



< e-^^^ (6.6) 



10 



We take here daV = V. 
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with < /c < 5. For the magnetization m{h, L) and its derivative, the susceptibihty 
x{h, L) = dm{h, L)/dh, these bounds yield 

h) = ijM„(L, h) + ij^^li^ ta^(-^£^) + 0(e--) (6.7) 

and 

x(L,/.) = ^Xo(i,M + ^^ta^(-^) 

+ J,(^)%o.,-V^).0(e-). (6.8) 

Here xo{L,h) = dMo(L,h)/dh and Ax{L,h) = dAM{L,h)/dh . 

In order to obtain Theorem A, and more generally the corrections to it in terms of an 
asymptotic power series in we proceed in several steps: 

i) We expand the functions AF(L, h), Mo{L, h), AM{L, h), Xo{L, h) and Ax(L, h) around 
the point ht{L) where AF{L,h) = 0, obtaining a power series in {h — ht{L)) with 
coefficients that are derivatives of AF{L, h) and Fq{L, h) at the point ht{L). 

ii) We Taylor expand the coefficients in i) into a power series in {ht{L) — ht). Combined 
with the volume, surface, corner expansion for the derivatives of F±{L,h) and the 
fact that ht{L) — ht can be represented as an asymptotic expansion in powers of 1/L, we 
obtain the coefficients of i) as power series in 1/L, with coefficients that are derivatives 
of the infinite volume free energies f±{h), surface free energies f^~^\h), and corner 
free energies f^\h) at the infinite volume transition point ht- 

iii) At ht, the derivatives of f±\h) are identified with the one-sided derivatives of the free 
energies defined by (2.15). 

iv) We use Lemma 6.1 below to replace the argument of the hyperbolic functions in (6.7) 
and (6.8) by few expansion terms with an additive error. 

v) In a final step, we use Lemma 6.3 to replace h — ht{L) hy h — h^{L), where h^{L) is 
the position of the susceptibility maximum. 

Lemma 6.1. Let x and y two nonzero real numbers which have the same sign. Then 

I 1 I . .tanhx tanhy ^ \ 

tanha; — tanhy < mm( , )\x — y\ (6.9a) 

X y 

and 

I 2 ,_2 I . .tanhx tanhy 

cosh a; — cosh j/ <2mm( , )\x — y\. (6.9b) 

' X y 

Lemma 6.2. For large L there exists a unique point ht{L) e U for which 
F^{L, h) = F_{L, h). This point satisfies the bound 
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Lemma 6.3. For large L there exists a unique point h^{L) & U as well as a unique 
point hu{L) e hi for which the susceptibility x{L,h) and the Binder cummulant U{L,h), 
respectively, attain its maximum. To the leading order in 1/L, their shift with respect of 
the point ht{L) is given by 



and 



h^{L) ^ ht{L) + 6 



hu{L) = ht{L) + 4 



^+-^- - + 0(—) 

(m+ -m_)3 L^d ^^^L^d) 



(m+-m_)3 L^d^ ^L^d)- 



(6.11) 



(6.12) 



Proof of Theorem A. Let us begin with the identification iii). Introducing 

d'fL'\h) 



my{h) = — — , x±{h) 



dh 



dK^ 



■, k = d,...,0, 



we get 



and 



f^^\ht) = hm /('^V/i) 



df^^\h) 



dh 



ht±0 



d''f^^\h) 



dh^ 



ht±0 



(6.13) 
(6.14) 
(6.15) 



for /c = (i, . . . , 0. In particular, the one-sided derivatives (2.16) as well as the limits (2.17) 
and (2.18) are expressed in terms of derivatives and limits of the corresponding differen- 

(k) 

tiable functions : 



and 
Also 



dh' 

,{d-i), 



h=ht 



r± = fr'ih). 



AF{L) = AF{L,ht) 



(6.16) 

(6.17) 
(6.18) 



with AF(L, h) = {F+{L, h) — F+(L, h)). To show all this, we first notice that, for the two 
phase case, the bound (3.19) from Theorem 3.1 reads 



d'' 



dh'' 



[Z{L,h) 



^-F+{L,h) _ ^-F_{L,h) 



< \V\''+^ max(e-^+(^''^), e-^-(^''*))0(e-^^^). 

(6.19) 

Taking into account that F±{L, h) are asymptotically dominated by f±{h)L^, the bound 
(6.19) implies that for h > ht the free energies f^''\h), k = d,...,0, defined by (2.15) 
actually equal the corresponding free energies f^\h), k = d, . . . ,0, from Theorem 3.1 (we 
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have chosen the notation for which min(/_|_(/i), /_(/i)) = f+{h) for h > hf). Similarly, 

f^'^\h) = f^\h)j k = d, . . . ,0, for every h < hf.^^ This identification immediately implies 
the equalities (6.14)-(6.18). 

Notice also that by (3.9) and Theorem 3.1 iii) one has 



dh'' 



< Co + e 



-6f 



(6.20) 



and thus also 



m±(/it(L)) = m± + 0{\ht{L) - ht\) ^ m± + O(i^tM) 



(6.21) 



according to Lemma 6.2, where we evaluate AF{L) with the help of (2.8) and Theorem 
3.1 iii). 

Expanding now M±{L, h) and F±{L, h), we have 



M±(L, h) = M±{L, ht{L)) + 0{{h - htiL))L^) 

= m±{ht{L))L^ + 0{{h - ht{L))L'') + 0{L''-^) 
= m±L'^ + 0{{h - ht{L))L'^) + 0((1 + \\n\\)L'^-^) 



(6.22) 



and 



-AF{L, h) = (^M+(L, ht{L)) M_iL, ht{L)^ (h - ht{L)) + 0{{h - ht{L)fL'') 

= (m+ - m_){h - htiL))L^ + 0{{h - ht{L)fL^) + 0((1 + ||«||)(/i - ht{L))L''-^) 
= 2x{l + 0{{h - ht{L))) + 0((1 + ||«||)L-i) . (6.23) 



Here 



(6.24) 



Using Lemma 6.1 to replace the argument of the hyperbolic functions in (6.7) and (6.8) 

by x, we get 

m(L, h) = '^+ + '^- + ^+~^- tanh(a;) + 0{{1 + \\k\\)L-') + 0{{h - ht{L))) (6.25) 



and 

X{L, h) 



m.+ — iri- 



cosh-"(x)L^ + 0((1 + \\k\\)L'^-^) + 0{{h - ht{L))L'') . (6.26) 



■"^^For h = ht, the asymptotic behavior will be determined by the first k = d — 1, ... ,0, for which 
fi^\ht) 7^ £^\ht). For example, if /^'-'^(/it) > fi^-^\ht), then /('=)(/it) = fi^\ht) for all k = d, . . . ,0 
{ofcomse, ff{ht) = f^^\ht)). 
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In order to replace further the argument x of the hyperbohc functions by the argument 



X = 



ih-h^iL)W 



used in (2.19) and (2.20), we finally use Lemma 6.3 to bound 

I tanhx — tanh£| < |a; — £| < 0{L~'^) 



and 



cosh ^ a; — cosh ^ x\ < \x — x\ < 0{L '*) 



(6.27) 

(6.28a) 
(6.28b) 



Combining the bounds (6.24) and (6.28) with the assumption \h — h^{L)\ < 0((1 + 
||k||)L-1), we get the bounds (2.19) and (2.20). 

The shifts (2.21a) and (2.21b) as well as the bound (2.22) on the mutual shift hu{L) — 
h^{L) follow from Lemma 6.2 and 6.3. □ 

Proof of Lemma 6.1. Without loss of generality, we may assume that a; > y > 0. Since 
is a decreasing function of t, we have > -^^^ and thus 

t * ^ y X 



tanh X — tanh y 



tanha; 



tanh y y 
1 — ^ < 1 - - 



tanh a; 



X 



\x - y\ 
\x\ 



(6.29) 



This concludes the proof of (6.9a). In order to prove (6.9b), we bound 



cosh ^ X — cosh ^ y 



-2 





fy sinht , 




r 2 


= 2 


/ 3-^^^ 


< 2 


/ cosh t dt 




Jx cosh t 




J X 



2 |tanha; — tanh j/| 

(6.30) 



and use (6.9a). □ 

Proof of Lemma 6.2. Using the bounds (6.2) we get, for sufficiently large L, the bound 



\ - dh - 



(6.31) 



Since AF(L) = AF{L,ht) = 0{L'^-^), we get the existence of a unique ht{L) for which 
AF(L, h) = 0. Moreover, ht{L) e {h -^.ht + f ) for some B. For h in this interval, the 
Taylor expansion of AF{L, h) around ht yields 

AF(L, h) = AF(L) - (m+ - m-)L'^{h - ht) + (h - ht)0{L'^-'^) (6.32) 

whenever L > B. This implies (6.7) (valid also for AF{L) = when ht{L) = ht). □ 

Proof of Lemma 6.3. To get (6.11), we may actually follow the proof of Theorem (3.3) in 
[BK90], replacing only ht by ht{L). Thus we use first (6.8) combined with the bound 



d^F±{L,h) 
dp 



<CkL'' 



(6.33) 
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(here may be chosen Ck = (Co)'' + 0(e-^^) by (3.9) and Theorem 3.1 hi)) to get 



\xiL,h) - ^{AM{L,h)fcosh-\^^^^)\ < 1 + 2(72. 



(6.34) 



Hence 



X{L, ht{L)) - x{L, h)>^ (AM(L, ht{L)))'- 



^(AM(L, h)Y cosh-^(^^^) - 2 - 4C,. (6.35) 



Next, we use (6.31) and (6.33) to bound \-^{AM{L,h))'^\ by S^CsL^^. As a consequence, 
\{AM{L,h)f -{AM{L,ht{L))f\ < 8AC2L^'^\h - ht{L)\ (6.36) 



and 



X(L, h,{L)) - x(L, h)>^ (AM(L, h,{L)))\l - cosh'^ (^^^)) - 



2AC2\h - ht{L)\L'' cosh-^{^^^%^ 



) - 2 - 4C2. (6.37) 



On the other hand. 



a 



h-ht{L)\L'^ < 



AF{L,h) 



< A\h - ht{L)\L' 



(6.38) 



by (6.31) and the fact that AF{L,ht{L)) = 0. Using the lower bound 

2 

cosh^ a> (1 + ^)^ > 1 + CK^ > 2\a\ 
vahd for any a, we imply that 

_2,AF(L,/i). 



(6.39) 



cosh (- 



-)< 



2 ' - a\h-ht{L)\L'^' 
Thus, using once more (6.31) and (6.38), we get 



(6.40) 



4Ld4" ■ 



x(L,/.,(L))-x(L,/i)>^(-L^)'(l-cosh-2(-i?)) 



a 



2 - 4C2 (6.41) 



whenever we suppose that \h—ht{V)\ > BL ^, where 5 > will be chosen later. Observing 
that 



cosh"^(-S) < 1, 
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it is clear that the right hand side of (6.41) is positive, once L is sufficiently large (how 
large depends on the choice of B). 

Thus, it remains to consider the case \h — ht{L)\ < BL~'^. Taking into account that 
AF{L,ht{L)) = we get 



dxiL,h) 



dh 



1 d^Mo{L,h) 



ht{L) 



dK^ 
+ 



+ 



ht{L) 

3 dAM{L,h) 



dh 



AM{L,ht{L)) + 0{e-^^^). (6.43) 



ht{L) 



Using the bound (6.33) we get 
dx{L,h), 



(6.44) 



Applying once more the bound (6.6), this time for k = 3, and using (6.31) we get, for 
\h-ht{L)\< BL-d, the bound 



d^yiL.h) ,1 , , . 1 — 3 tanh^(-^^^7r^) . 



(6.45) 



Taking into account that according to (6.38) one has |AF(L, h)\ < 2A\h — ht{L)\Ld and 
choosing i? > so that 

1 - 3taiili-(_4Z?) 



£ := 



cosh^(^5) 



>0, 



we get 



(6.46) 



(6.47) 



for \h — ht{L)\ < BL ^ and L large enough. The bound (6.47) together with the bound 
(6.44) implies that there exists a unique zero h^{L) of in the interval {ht{L) - 

-gr, h{L) + ^) and that h^{L) - ht{L) = ©(L'^'^). 

For h - ht{L) = 0(-^) we have AF{L,h) = 0{^) by (6.38) and thus, using (6.45) 
and (6.36) we get 



- -2(^(m+ - m_))^L3'^ + 0(1^^-'). 



(6.48) 



Taking into account (6.44) we conclude the bound (6.11). 

To prove the bound (6.12), we first notice, by a straightforward computation, that 



(6.49) 
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Using the fact that, for L large, 

(Px_ 1 d^Mo{L,h) 1 d^AM{L,h) AF{L,h) 

dh? ~ ~1? rfP ^ 21? dh? ^ 2 

1 d^AM{L,h) AM(L,/i) 3 , dAM{L,h) .2 

Ard\ IT. / 



dh^ cosh2(^^) 4L<^' d/i ^ cosh2(^^) 



<Q(i^-)-i^-^(g)- ^-": iL;7i-' . (6.50) 



we find that U{L,h) is negative (and thus smaller than U{ht{L)) whenever 

l_3tanh^(^flii)) 



for some positive £. To meet this condition, it suffices to take \h — ht{L)\ > with B 
such that 



2,aB, _ 3 



cosh^(— )>-(! + 6) (6.52) 

for some e > 0. Indeed, using (6.38) we get |AF(L, /i)| > ^L'^\h — ht{L)\ and thus 
cosh2(^^^^) > 1(1 + e), which implies (6.51) with e = 2e/(l + e). 

In the interval \h — ht{L)\ < ^ we consider the leading terms to ^^^^'^^ and ^Jh^'^^ • 
Namely, 



dh 3 ^ 2 2 ^ 2 ^ 

dAM{L,h) \Tp/T U\ dAMo{L,h) 



dh 



,„.,3,A£lM),_3=i=.n,,^))] (0.53) 



yielding 



AM{L,h) ' 2 ' AM{L,h) 
dU{L,h) 



dh 

and 

d2C/(L,/^) 2 _ ,2,^ , ^ .,AF(L,/.) 



h,{L) 3 m+-m_ 



--L^^(m+ -m_)2(l + 2sinh^(^=^^^f^)) + 0(L3'^-i) < 



d/t2 3 



< - jL^'^(m+ - m-f. (6.55) 
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Thus, there exists a unique root hu{L) of the equation dU{L,h)/dh = in the interval 
\h - ht{L)\ < ^ and hu{L) - ht{L) = 0{L-'^'^). Moreover, for h - ht{L) = 0{L-^'^) 
we have AF(L, h) = 0{L-<^) and thus (PU{L,h)/dh'^ = -|L3rf(m+ - m_)2 + OiL^"^-^). 
Hence, 



dU{L, h) dU{L, h) 



dh dh 



d'^U{L, h) 



+ 



\h-ht{L)\ = Q (6.56) 



ht{L)+S,{h-ht{L)) 

yields the shift (6.12) claimed in the lemma. □ 
We are left with the proof of (2.23) and (2.24) for 



|ft_ft^(i)|>i^!(LdMI. (6.57) 

' ^' (m+ - m-)L ^ ' 



Using (6.2), (6.11) and the bound (6.38), we first note that the condition (6.57) implies 

(6.58) 

Combined with (6.7) and (6.8) we obtain 

m(L, h) = L-'^M+{L, h) + 0(e-"^/^) 
X(L, h) = L-'^^M+^lR + 0(e-"^/^) 

m(L, h) = L-'^M_{L, h) + 0(e-"^/^) 

ii h < h^{L) — (^^Tz^ii^^- Expanding M±(L,h) and its derivative into volume, surface, 
corner terms, this leads to 

m{L,h) = m+{h) +0{1/L) 
Xmh)^X+{h) + 0{l/L) 



if^>M^)+(^™ji:and 



m{L,h) =m-{h) + 0{l/L) 
x{L,h)=x-{h) + 0{l/L) 
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if/i</i.(L) + 7Mi«I 
aV / {m^—m_)L 

Next, we recall |e±(c) — e±\ < \\k\\ for the asymmetric Ising model (2.1). As a con- 
sequence, |AF(L)| < 2dL'^-'^{2\\K\ \ + 0{e-^^)). Combined with Lemmas 6.2 and 6.3 we 
conclude that 

IM^-fe.|< ,"""''"^;' (6.58) 

(|m+ — m-\)L 

As a consequence, h < h^{L) + ^^^"'"J^^^^j^ implies h < ht and hence m{h) = m-{h) and 

X{h) = while h > h^{L) + implies h > ht and hence m{h) = mj^{h) and 

x{h) = x+(/i). The condition (6.57) therefore implies the bounds (2.23) and (2.24). 
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Appendix A: Strong isoperimetric inequality 
Using the standard isoperimetric inequality, 

\dW\> i cil^l^, (A.l) 

in the proof of Lemma B.3 below, we would get, for (i > 4, a negative factor on the 
right hand side of the bound (B.2). We strengthen (A.l) with the help of an additional 
information — the fact that considered sets W are finite unions of closed elementary cubes. 

Lemma A.l. Let W he a union of closed elementary cubes. Then 

\dW\>2d\W\^ . (A.2) 



Proof. The proof is just a particularly simple case of the proof of optimality of the Wulff 
shape [T87]. Namely, 

\W + eC\ - \W\ 
\dW\ = lim ' ^ ' (A.3) 

e^O e 

where eC is the rescaling, by the factor £, of the (hyper)cube C of side 2 with the center 
at the origin, and 

W + eC = {x + y : X eW,y e eC} (A.4) 

is the ^-neighborhood of W in the maximum metric. The Br unn- Minkowski inequality 
(valid also for nonconvex W see for example [Fe69]) yields 

\W + eC\^ >\W\^ + \eC\^. (A.5) 

Thus 



\WA-eC\-\W\ i\W\^ + \eC\^Y - \W\ i 
lim ' ^ ' ^ > lim ^' ' ' ' ' ^ = d\C\-^\W\— = 2d\W\ — . (A.6) 



d-l , , d~l 



□ 
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Appendix B: Proof of Lemmas 4.1 - 4.5 



We start with two Lemmas, B.l and B.2, that are an important technical ingredient to 
prove Lemmas 4.1, 4.2 and 4.4. 

Lemma B.l. Let Y be a short contour with suppF Pi dV C dK{k). Then 

i) supp YndV = supp Y n dK{k) = supp Yn{dVn dK{k)) 

ii) a Int y C dK{k) U d supp Y 

Hi) Int Y ndV = IntY n dK{k) = Int y n {dV n dK{k)) 

Proof, i) suppy n dK{k) = suppy n dK{k) nV = suppy n dK{k) n dV since V n 
dK{k) = dV n dK{k) and supp y C On the other hand supp Y fl 91/ C dK{k) imphes 
supp y n C dK{k) n 9y and hence supp Y fl 91/ C supp Y n dK{k) n Combining 
this with supp Y n dK{k) fl 9F C supp y fl dV we obtain i). 

ii) Follows from the fact that all components of Int Y are components of K{k) \ supp Y . 

iii) In order to prove iii), we first prove Inty fl dV C 9Inty fl dV . This can be 
proven as follows: the inclusion Inty C V implies C (Inty)"^ implies dist(x, > 
dist(a;, (Int y)'^). Therefore dist(x, (Inty)^) = for all x G dV and hence for all x G 
Inty n dV . Since x G Inty and dist(a;, (Inty)^) = implies x G 9 Inty, this proves 
Inty nay C d\niYr\dV. 

Using ii), the (just proven) fact that Inty fl dV C 9 Inty fl dV and the fact that 
asuppynSy C dK{k), one proves that IntY (IdV C dK{k). Intersecting both sides with 
Int y and observing that Int y C F while V n dK{k) — dK{k) n dV , one concludes that 



proves iii). □ 

Lemma B.2. Let Yi and Y2 be two non- overlapping contours with suppyi C Inty2. 

Assume that Yi is a short contour with dV ^ suppy2 C dK{k) some corner k. Then Yi is 
a short contour with dV fl supp Yi C dK{k) as well. In addition supp Y2 fl Int yi = 0. 

Proof. By Lemma B.l iii) and the fact that suppyi C Inty2, suppyi H dV C Inty2 fl 
dK{k) C dK{k). Let now xq G suppyi C Inty2. By the definition of Inty2 and Exty2, 
and by the fact that suppy2 is connected, we may construct a path ui in K{k) which 
connects xq = a;(0) to infinity such that 



Int Y ndv GintY n dK{k) = Int y n dK{k) n dV. 



This combined with the fact that 



Int y n dK{k) n c Int y n av. 



uj{t) G Inty2 
Uj{t) G suppy2 
oj{t) G Exty2 
uj{t) G K{k) \ V 



for t G [0, 1) 
for t G [1,2) 
for t G [2, 3) 
for t G [3, 00) 
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Assume now that supp Y2 fl Int Fi ^ 0. Since supp Y2 is a connected set which does not 
intersect supp Yi, this imphes that supp 12 C Int Yi. As a consequence X2 = a;(2) e IntYi 
and 0^1 [2,00) is a path in K{k) which connects X2 G Int Yi to infinity. But this imphes that 
a;|[2,cx)) must intersect the set supp Yi, and hence, by the assumption that supp Yi C Int Y2, 
the set IntY2. This is a contradiction, because u> was constructed in such a way that 
a;(t) ^IntY2 for t > 1. □ 

Proof of Lemma 4.1. 

i) Since supp Y2 C Ext Yi, supp Y2 flint Yi = 0. It foUows that each point in Int Yi can 
be connected to 5Int Yi (and therefore to supp Yi) by a path uj which does not intersect 
supp Y2. As a consequence, all points in Int Yi lay in the same connectivity component of 
V \ supp Y2 as supp Yi. Since supp Yi C Ext Y2, we concluded that Int Yi C Ext Y2, and 
hence Int Yi U supp Yi C Ext Y2. In a similar way, Int Y2 U supp Y2 C Ext Yi. 

ii) Let us first assume that Y2 is a short contour. Then Yi is a small contour as well 
and supp Y2 fl Int Yi = by Lemma B.2. Continuing as in the proof of i) we obtain that 
Int Yi UsuppYi C C2. 

Next, consider the case where Yi is short while Y2 is long, and assume that supp Y2 fl 
Int Yi ^ 0. Since supp Y2 is connected, this would imply that supp Y2 C Int Yi. By Lemma 
B.2, this would imply that Y2 is short as well. Therefore supp Y2 fl Int Yi must be empty. 
Again, this implies Int Yi U supp Yi C C2. 

As the last case, assume now that both Yi and Y2 are long. Since C2 C Int Y2, IC2I < 
by the definition of the exterior for long contours. Since supp Y2 is a connected set, 
while C2 is a connected component of F \ supp Y2 , both C2 and V \ C2 are connected 
sets. Observing that supp Yi C C2 implies V \ supp Yi D F \ C2, we then introduce the 
component Ci of V \ supp Yi which contains V \ C2. A minute of reflection shows that 
|Ci| > IV" \ C2I, which, by the fact that |V" \ C2I > implies that Ci = Ext Yi. As a 

consequence 

IntYi C F\ExtYi C C2, 

which concludes the proof of ii). 

iii) follows from ii). □ 

Proof of Lemma 4.2. 

We only have to show that supp Yi C Int Y2 and supp Y2 C Int Yi leads to a contradic- 
tion. In fact, suppYi C IntY2 implies that suppYi UlntYi C IntY2 by Lemma 4.1. As 
a consequence Ext Yi D supp Y2 U Int Y2 which implies that supp Y2 C Ext Yi . But this is 
incompatible with supp Y2 C Int Yi . □ 

Proof of Lemma 4.4. 

Let Yfc be an internal contour. Due to Lemma 4.2, this implies that (Int Y^ U supp Y^) C 
Int Yj for some j ^ k and hence Ext = \ Ui^^fcl-^^^ Usupp Yi). Iterating this argument, 
we get that the set Ext is given by 



n 

Ext = F \ y (Int Y[ U supp Y.^), 



(B.l) 
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where {Yf,...,Y^} are the external contours in {Yi, . . . , Obviously, Ext is separated 
from the rest of V by the support of the contours Yf,...,Y^. We therefore only have to 
show that Ext is connected. 

Let Ei = V\ (Int Yf U supp Yf) = Ext Yf and Ek = Ek-i\ (Int Y,^ U supp Y,^). Assume 
by induction that Ek-i is connected. Let x,y & Ek- We have to show that x and y can be 
connected by a path cok in Ek- Using the inductive assumption, x and y can be connected 
by a path oJk-i in Ek-i- Assume without loss of generality that oJk-i intersects the set 
W = Int Yj^f U supp y^f, and let xi be the first and yi the last intersection point of u>k-i 
with W. Since both W and V \ W are connected, the boundary 

dvW = {xeV\ dist(a;, W) = dist{x, V\W) = {)} 

is connected, and xi and yi can be connected by a path oj in dyW. Using the path Uk-i 
from X to xi, the path lo from xi to yi, and again the path Uk-i from yi to y, we obtain a 
path (Dfc in E^ U 9y which connects x to y. The desired path ojk is obtained by a small 
deformation of uik which ensures that is a path in E'^. □ 

In order to prove Lemma 4.5, we need the following Lemma, which is based on the 
strong isoperimetric inequality proven in Appendix A. 

Lemma B.3. Let W be a union of elementary cubes in V , with \W\ < Then 

\dW n dV\ < lud^] \^^ \ (^-2) 

^ ' J. 

\dW\ <l^l + ^17^ J \dW \ dV\ (B.3) 

Proof. We introduce the {d — 1) dimensional faces 

Fi = {x eE."^ \xi = 1/2, 1/2 <Xi<L + 1/2} i = l,... ,d 

Fd+i = {x eR"^ \xi = L + 1/2, 1/2 <Xi < L + 1/2} i = l,... ,d, 

together with the projections ni : V —>■ Fi, iTd+i ■ V — > Fd+i, where x' — 7ri{x) has 
coordinates x'f^ = Xk ior k ^ i and x'^ = 1/2, while x'^ = L + 1/2 for x' = 7rd+i{x). Finally, 
for each elementary {d — l)-cell p e C, we define tt-(p) as the projection 7rj(p) onto the 
face Fi which is parallel to p, and 7r_|_(p) as TTd+iip)- 

Let Gi = Fi n dW , i = 1, . . . ,2d, and consider an elementary {d — l)-cell p G Gi, 
together with the line £ that links the center of 7r_ {p) to the center of n+ {p) . Then £ must 
intersect dW an even number of times. Define 

Hi = {p ^ Gi \ there does not exist p' e dW \ dV with 7r_(p') = 7r_(p)} 

and consider an elementary {d — l)-cell p e Gi\Hi, i = 1, . . . ,2d. Then either both 7r_ (p) 
and 7r+{p) lay in U^liGj \ Hj, in which case there are at least two elementary {d — l)-cell 
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•p' e dW\dV with 7r_(p') = T^-{v)-, or only one of 7r_(p) and 7r_|_(p) lies in \J^^^Gj \ Hj, in 
which case there is at least one elementary {d — l)-cell p' e dW \ dV with 7r_ {p') = 7r_ (p). 
As a consequence, 

2d 

^\Gi\Hi\<\dW\dV\. (B.4) 
1=1 

On the other hand, 



by the fact that < Using the strong isoperimetric inequality (see Appendix A), 

we obtain that 

\Hi\ < ^2-^/'^\dW\. (B.5) 
Combining (B.4) and (B.5), we get 

2d 2d 2d 



i=l i=l i=l 

< \dW\dV\ + 2-^/'^\dW\ 

= (1 + 2-^''^)\dW \ dV\ + 2-^/'^|aW^ n dV\ 

and hence 

1 + 2"^/'^ 
\dWndV\ < — — -j:;\dW\dV\ 

which implies (B.2). The bound (B.3) follows from (B.2). □ 
Proof of Lemma 4-5. We start with the observation that 

\Y\ = \Y\a+\Y\a-i + --- + \Y\i, (B.6a) 

where \Y\k denotes the number of A;-dimensional elementary cells in-^^ Y, and similarly for 
the boundary dWi of a component Wi of Int Y, 

\dW^\ = \dW^\d-i + ■■■ + \dW^\i. (B.6b) 

Using the fact that \dWi fl dV\d-i = \dWi n dV\d-i we then decompose \dWi\ as 

d-2 

\dWi\ = \dWi \ dV\d-i + \dWi n dV\d-i + J2 ■ (B.7) 

fe=l 



■"^^ As in Section 3, a A;-dimensional cell c in supp Y is only counted if there is no (fc + l)-dimensional 
cell c' in supp Y with c C c' . 
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For a long contour Y, we use Lemma B.3 applied to the set Wi, together with the fact 
that dWi C dWi to bound 

\dWi n dvU-1 < [^rjT^) \9Wi \ dvU-i. < [^Tjd^J 1^^^ \ ^vU-i- (b.8) 

For short contours Y, on the other hand, \dWi fl dV\ < \dWi \ dV\ which implies (B.8) 
with a better constant. Therefore (B.8) is valid for both long and short contours. 

As a last step we observe that each cube c in suppy can be shared by at most 2d 
elementary faces in {dWi \ V) U ■ ■ ■ U {dWn \ V), while each d—1 dimensional elementary 
face in Y may be shared by the boundary of at most two different components of Int Y U 
Ext Y. Since all lower dimensional elementary cells in Y belong to a unique component of 
Int y U Ext y , we get that 

n n d—2 

J2 \dWi \ dV\d-i + l^^^l'^ ^ 2d|y|. (B.9) 

i=l i=l fc=l 

Combining (B.7) with (B.8) and (B.9) and the fact that 

n 

Nqv {Int Y) < J2 \9Wi ndV\, (B.IO) 

we obtain the first two bounds of the lemma. 

In order to prove (4.4) we observe that V{Y) — suppy U Wi U ■ • • U Wn which in turn 
implies dV{Y) C d supp Y U dWi U • • • U dWn and hence 

\dV(Y)\ < |asuppy| + \dWi\ + ■■■ + \dWn\ 

Combined with the bound |9suppy| < 2d\Y\ we obtain the remaining bound of Lemma 
4.5. □ 
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Appendix C: Inductive Proof of Lemma 4.6 



In this appendix, we prove Lemma 4.6. Actually, we will first prove the following Lemma 
1. In order to state the lemma, we re 
auxiliary contour model with activities 



C.l. In order to state the lemma, we recall the definition of fq^^ as the free energy of an 



, . , r K'(Yi) if \V(Yi)\ <n, , , 

[ otherwise, 

and define 

/(")=min/(-), (C.2) 

4"^ = - /• (C.3) 

We also assign a number v{W) to each volume of the form (4.7a), 

v{W) = m■^.^\V{Y)l (C.4) 

YmW 

where the maximum goes over all contours Y with suppF C W. Obviously, t'(Inty) < 
for all contours Y. In fact, 

v{IiitY) <\V{Y)\ (C.5) 

due to the fact that dist(y , y) > 1 if y is a contour in Int Y. 

Finally, we recall that for a volume W of the form (4.7a), \W\ is used to denote the 
euclidean volume of W, while for a contour Y and the boundary dW of a volume W, \Y\ 
and \dW\ are used to denote the number of elementary cells in Y and dW, respectively 
(see Equation (B.6a) and (B.6b)). 

Lemma C.l. Let 

, = e-r(i-(2C,+i)7)ga+2 - ^ ia-2)2d^ 

Then there is a constant eo, depending only on d and N, such that the following statements 
are true for all e < eo and all n > 0, provided \V{Y)\ < n, v{W) < n, and o: > 1. 

^) \K'^{Y)\ < e\^\ . 
\) If a^q^\V{Y)\^''^ < a then ^JY) = 1. 



iii) If aq-^\V{Y)\^/'^ <a then K'^(Y) = Kq{Y) . 

iv) If aq^\W\^/'^ <a then Zq{W,h) ^ Z'^{W,h). 



Proof. We proceed by induction on n, first proving the lemma for n = and then for any 
given n e N, assuming that it has been already proven for all integers smaller than n. 
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Proof of Lemma C.l for n=0. 

For \V{Y)\ = we have XqiY) = 1 and thus K'g{Y) = Kq{Y) = p{Y). This makes i), 
ii) and iii) trivial statements. Using iii) for = 0, we then conclude that Zg{W, h) = 

Z'g{W, h) for v{W) = 0. By i) Z^^W, h) = Z'q{W, h) and thus the partition function can be 
analyzed by a convergent expansion yielding 

Observing that > f^^\ this concludes the proof of Lemma C.l for n = 0. 
Proof of Lemma C.l i) for = n. 

Due to (C.5), v(Inty) < n, and all contours Y contributing to Z^(Int^ y, /i) obey the 

condition |F(y)| < n. This implies that |Kq(y)| < e'^' by the inductive assumption i). As 
a consequence, the logarithm of Z^(Int^y, h) can be analyzed by a convergent expansion, 
and 

logZ;(Int^y,/i) +/("-i)|Int^y|| < 0(e)|aint^y| +7TiVav(Int^y). (C.7) 
Combining (C.7) with the induction assumption v), we get 



n 



Z-m 

(Int^y, h) 



Z'JlntmY,h) 



< ga("-i)|Inty|g(2Ci7r+0(e))|y|^ (^(j g^ 

where we have used Lemma 4.5 in the last step. Observing that 

\em~f^-'^\<0{e), (C.9a) 

which implies the bound 

|(e,-eo)-4"-i)| <0(e), (C.9b) 
we use the assumptions (3.7) and (3.11) to bound 

|p(y)e^«(^)| < e-^l^le^^^s^('"PP^)e(^«-^°)l^l'' < e-^^-^^-^^'^^l^le''^""'^!^!'*. (C.IO) 

Here |y|d is defined as the number of d-cells in Y and thus |V(y)| = |Inty| + |y Id- 
Combining now (C.IO) with (C.8), we obtain 

\K'(Y)\ < (y)e4"-^^l^mie-(--0(^)-(i+2^^)^-)l^l. (C.ll) 
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Without loss of generahty, we may now assume that Xg(^) > (otherwise K'^(Y) = and 
the statement i) is trivial). By the definition of x'qiY), this implies 

(/f - f^-'^)\ViY)\ < 1 + < (1 + a)\Y\ 
for all m 7^ g. As a consequence, 

4"-i)|F(y)|<(i + a)|y|, (C.12) 

provided Xq{Y) 7^ 0. Combined with (C.ll) and the fact that Xq{Y) < 1, this implies that 

\mY)\ < e-[--l-O(e)-a-(l+2C07r]|r|^ (^.13) 

which yields the desired bound i) for |F(y)| = n. 

Proof of Lemma C.l ii) for k = \V{Y)\ < n and a'^''^\V{Y)\^/'^ < a . 

We just have proved that i) is true for all contours Y with |F(y)| < n. As a consequence, 
both fm^ and fm^ may be analyzed by a convergent cluster expansion. On the other hand, 

by the isoperimetric inequality and Lemma 4.5. Using this bound and the definition of 
fm \ one may easily see that all contours Y contributing to the cluster expansion of the 
difference fm^ — fm^ obey the bound 

|y|>^(/, + l)(d-i)/d>|ffci/d=: no. 

As a consequence, 

\fi^-f^^\<{Ker'', 

where X is a constant depending only on the dimension d and the number of phases N. 
Using the bound (C.14) for the second time and recalling that |V^(y)| = k, we get 

\fi^ - f^^mY)\ < (i^e)""|y(y)|V^g|y| = 0{l)no{Ker>\Y\ < 0(e)|y|. (C.15) 

Combining (C.15) with the assumption a^J^^\V{Y)\^/^ < a and the bound (C.14), we obtain 
the lower bound 

o^\Y\ - [/f ) - f^^]\ViY)\ > a\Y\ - a(-)|F(y)| - Oie)\Y\ > 

0(e)) \Y\ = (2 - 0{e))\Y\ > 2 - 0(e), 

where, in the next to the last step, we used the definition of a, see (C.6). Combined with 
(4.16b) we obtain the equality x'qiY) = 1. 
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Proof of Lemma C.l iii) and iv). 

Given ii) and the definition of K'^iY) and Z'^{W,h), the statement is obvious. See 
[BK90] for a formal proof using induction on the subvolumes of W and Int Y. 

Proof of Lemma C.l v). 

We say a contours Y is small if a^J''^\V{Y)\^/'^ < a while it is large if a^J''^\V{Y)\^/'^ > a. 
We then use the relation (4.8) to rewrite Zq{W,h) by splitting the set {Yi,--- ,YA;}ext 
of external contours into {Xi, ■ ■ ■ , Xf} U {Zi, • • ■ , Z^"}, where Zi, - ■ ■ , Zk" arc the small 
contours in {Yi, ■ ■ ■ , Yfcjext and Xi, ■ ■ ■ , Xk' are the large contours in {Yi, ■ ■ ■ , Yfcjext- Note 
that for a fixed set {Xi, ■ ■ ■ Xk'}, the sum over {Zi, ■ ■ ■ Z^"} runs over sets of mutually 

k' 

external small q'-contours in Ext = U V{Xi). Resumming the small contours, we thus 

i=l 

obtain 



k' 



Zg{W, h)= YI ^r""(Ext, h) H p{Xi) \[ ZmilntmXi, h) 

{Xi,...,Xfc,}ext i=l 



(C.16) 



Here the sum goes over sets of mutually external large contours in W and Zg™^^^(Ext, K) 
is obtained from Zq(Ext, h) by dropping all large external qf-contours. 

Due to the inductive assumption iii), Kq{Y) = Kq{Y) if y is small. Since [^^(y)! < e'^' 
by i), Z^™*^'(Ext, h) can be controlled by a convergent cluster expansion, and 

|Z^-^"(Ext,/i)| < e-/r"'|Ext|eO(e)|aExt|g7riVav(Ext)^ (C.17) 
where /g™^'^ is the free energy of the contour model with activities 

.small,^^ _ / KiY) if l^(^) I < ^ and y is smaU, 



^ ^ ^ 1 otherwise. ^ ^ 

On the other hand, 

ll\Zm{lntmX,,h)\ < e-/*"-^)|IntX,|gO(e)|aintX,|g7riV,,(IntX,) 
m 

by the induction assumption v). Observing that the smallest contours contributing to the 
difference of fm^ and fm~^^ obey the bound 

C3 63 
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while < n, we may continue as in the proof of (C.15) to bound 

l^(n-i) _ intX.I < 1/(^-1) - < 0{l)n^{Ker < 0{e) . 

Since | IntXj| < 0(l)|Xj| by Lemma 4.5, we conclude that 

ni^m(IntX,,/.)| <e-^'"'l^^*^^le^(^)l^^le^^^«-(^^^^^). (C.19) 

m 

Combining (C.17) and (C.19) with the bounds 

\p{Xi)\ < e-r\Xi\-eo\Xi\d^iTNsvi8uppXi) ^ g-(r-0(e))|Xi | g-/("> |X< |dg7TArev (supp X^) ^Q 20) 

and 

k' k' 

|aExt| < \dW\+^\dV{Xi)\ < \dW\+CsJ2\^i\^ (C-21) 
i=i i=i 

and the equality Ndv(E^i) + ^i=i[Ndv{su.pp Xi) + Ngvilnt Xi)] — Nqv{W), we conclude 
that 

k' 

X J2 e-[^r^"-^^"^ll^^*lf[e-("-0(^))l^^l. (C.22) 

{Xl,-..,Xfc,}ext »=1 

Next, we bound the difference /|™^^^ — fq^^. In a first step, we use the isoperimetric 
inequality together with Lemma 4.5 and the definition of large contours to bound 

1^1 > 7^19^(^)1 > ^iViX)]"^ > ^\V{Xr'^ > io := (^-23) 
Cs Cs 63 63 

for all large contours X. Next, we observe that 

|/f)-/r"l<(ire)'o<-^-j±^, (C.24) 
where K is a constant depending only on d and N. Recalling the condition a > 1, we get 

|^(n)_^small|< l^(n)^ (C.25) 

provided e is chosen small enough. Combining (C.22) with (C.25), we finally obtain 

|^g(W^,/i)| <eO(^)l^^le^^^«-(^)e-^'"'l^l J] g-^^l^^^' fj e"^!^^! (C.26) 
with 

f=(r-l). (C.27) 

At this point we need the following Lemma C.2, which is a variant of a lemma first proven 
in [Z84], see also [BI89]. 
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Lemma C.2. Consider an arbitrary contour functional Kq{Y) > 0, and let Zg be the 
partition function 

^'^W- E tl{K,{Y^e\^^^). (C.28) 

Let Sq be the corresponding free energy, and assume that Kq{Y) < e'^', where e is small 
(depending on N and d). Then for any a > —Sq the following bound is true 

J2 e«l^^*in^.(>^»)<e^^'^"''^', (C.29) 

{Yl,--- ,Yfe}ext » 

where the sum goes over sets of mutually external q-contours in W. 



In order to apply the lemma, we define Kq{Y) = e ^1^1 if F is a large q'-contour, and 
Kq{Y) — otherwise. With this choice, 

where £o is the constant from (C.23). As a consequence, 

(n) 

-S, < a := ^ (C.31) 

provided e is small enough. Applying Lemma C.2 to the right hand side of (C.26), and 
observing that e := e~'^ < e, we finally obtain the desired inequality 

This concludes the inductive proof of Lemma C.l. □ 

Proof of Lemma C.2. 

The partition function Zq is defined in terms of the polymer model with activities 
K*{Y) — ^q(y)el^'. For i small enough, Zq can be controlled by a convergent cluster 
expansion and 

llogZq(Inty) + Sq|Inty|| <0(e)|ainty| < 0(e)\Y\. 
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On the other hand, | W| = |Ext| + (| Int + 11"^!) if {Yi, • • • , Ifcjext is a set of mutually 
external contours in W. Combined with the fact that —a < Sq = 0{e), we obtain 

k 

< e'.m j2 nK,(r,)z,(intr,)eO(^')i^^i-^'''i^^i<^ 

k 

{Yl,---,Yfc}ext »=1 

Proof of Lemma 4.6. 

Lemma 4.6 i) through iv) follows from Lemma C.l and the fact that / = lim and 



ttg = lim ttq 

n— »oo 



in) 



In order to prove the statement v) , we extract the factor 

(") (") 
max e-^|E^*le-5^il^^l < max e"^! 'e-^^/^^^) E ja^x.)! 

{Xi,---,X^/} {Xi,--- ,Xf,,} 

(n) 

< max e-Vl^\t^le-(-/4C73)|a[/| 

UCW 

from the right hand side of (C.26), and bound the remaining sum as before. Taking the 
limit n — > oo in the resulting bound, this yields 

\Zq(W,h)\ < e^-^avWe(0(e)+0(e-3V*))|9W^|g-/|M^| ^\W\U\ ^-(r /^Cs)\dU\ ((..32) 

We conclude, with the help of the isopcrimetric inequality, that 

\Z (W h)\ < e^^^s^^^)e^^^')+^(^"'"^'^)l^^le--''l^l max e-^l^\^l-(2dT/4C3)|[/|'*/(<*-i) 

' ' ~ UCW 

where we used the fact that the maximum over U is obtained for either U = W or U = (l). 

Observing that 2d\V\'^/^'^-^^ = \dV\ and that Ndv{V) can be bounded by \dV\, the 
bound (C.33) implies Lemma 4.6 v). □ 

Proof of the bound (4.26). 

Due to the bound (C.12) we have a';P~^^\V(Y)\ < (1 + a) 1^1 if Xg(^) 0. Using 
the strategy which was used to prove (C.15), we replace a^" by a^, concluding that 
Xq{Y) ^ implies aq\V{Y)\ < (1 + 0(e) + Q!)|r|. □ 
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Appendix D: Proof of Lemma 4.7 

We start with a combinatoric Lemma that will be used throughout this appendix. 

Lemma D.l. Let ko be a positive integer and let G{h) be a function which satisfies the 
bounds 

d^ 

-G{h) 



dh^ 



for all multi-indices k with 1 < < ko and some A > 0. Then 

d^ 



dh^ 



< |A;|!Al'=le«W 



for all multi-indices k with 1 < \k\ < ko. 
Proof. Observing that 



d" 



where Hk{h) is a polynomial of degree in the derivatives of G, the the Lemma is 
immediately obtained by induction on □ 

Keeping the notation of Appendix C, we now prove the following Lemma, which contains 
statements i) through iii) of Lemma 4.7. 

Lemma D.2. There is a constant K < oo, depending only on N , d and the constants 
introduced in (3.8), (3.9), and (4-16), such that the following statements are true provided 

€ < €o, a > 1 and n > 0. 

i) For \V{Y)\ < n and Hq eU one has 



dhf" 



< {Ke)\ 



(D.l) 



h=ho 



provided 1 < \k\ < 6. 
a) For v{W) < n and ho &U one has 



d^ 



dh^ 



\og Z'(W,h) 



(D.2) 



h=ho 



provided 1 < < 6. 
iii) For v{W) < n and ho &U one has 



d^ 



< |jt|!(Co + 0(e))|»'|)l'le--'l»'le°WI''»'le^'''<»'("'' (D.3) 



h=ho 



provided 1 < \k\ < 6. 
Proof. As in the proof of Lemma C.l we proceed by induction on n. 
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Proof of Lemma D.2 for n=0. 

For \V{Y)\ = 0, K'^{Y) = Kq{Y) = p{Y), which makes i) a trivial statement. As a 
consequence, the left hand side of (D.2) can be analyzed by a convergent cluster expansion, 
leading immediately to the bound (D.2) for v{W) = 0. Bounding finally (Cq^' + 0(e))|VF| 
by {(Co + 0(e))|W^|}l^l and observing that Zg{W,h) = Z'g{W,h) if v{W) = 0, we obtain 
iii) with the help of Lemma D.l. 

Proof of Lemma D.2 i) for |V(Y)| = n. 

Using the assumptions (3.8) and (3.9) together with Lemma D.l, the bound (C.IO) can 
be easily generalized to derivatives, giving 



p{Y)e 



< \k\\ [2Co|y|]l'=l ^-{r-ir-0{e))\Y\^a,\YU . 



(D.4)) 



In a similar way, the bound (C.8) can be generalized to derivatives, using the inductive 
assumptions ii) and iii) together with Lemma D.l and Lemma 4.5. This gives 



d^ 
dhP 



(Int^ F, h) 
(Int^r,/i) 



< |A;|![(2Co + 0(e))|Inty|]l'^le"<'l^"*^le('^i^"+^(^))l^l. (D.5) 



Using finally the possibility to analyze the derivatives of fm ^\h) by a convergent expan- 
sion due to the inductive assumption i), we bound 



d'^f^^-'\h) 



dh^ 



<Ci''+0(e)<(Co + 0(e))l 



As a consequence. 



d^ 



dh^ 



X'{Y) 



< {C^\V{Y)\)\'^\ 



(D.6) 



for all multi-indices of order \k\ < 6. Here Ci is a constant that depends on N and 
the constants introduced in (3.8), (3.9) and (4.16). Combining (D.4) through (D.6) and 
bounding terms of the form 0(1)|F(F)| and 0(1)|F| by e'^^-*^)'^', we obtain the bound 
(D.l). 

Proof of Lemma D.2 ii) for v(W)=n. 

We just have proved that i) is true for all contours Y with |V^(1^)| < n. As a consequence 
the derivatives of log Z'q{W,h) can be analyzed by a convergent cluster expansion. The 
bound (D.2) immediately follows. 
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Proof of Lemma D.2 v) for v(W)=n. 



We define: a contours Y is small if aq{ho)\V{Y)\^/^ < a, while a contour Y is called 
large if aq{ho)\V{Y)\^/'^ > a. As in Appendix C we then rewrite Zq{W, h) as 



Zq{W, h)= J2 ^r""(Ext, h) n P{Xi) n Zm{ln.tmX,, h) 



(D.7) 



where the sum goes over sets of mutually external large contours in W and Z^™ (Ext, h) 
is obtained from Zq(Ext, h) by dropping all large external q'-contours. 

Due to Lemma 4.6, Kq(Y) = K'q{Y) if Y is small and h = ho. Combining this with the 
bound (4.13) and the inductive assumption (D.l), we conclude that 



dh 



< (2Ke)l 



(D.8) 



in a certain neighborhood Uq of ho. As a consequence, the derivatives of log ^^"'^'^(Ext, h) 
can be controlled by a convergent cluster expansion, and 



d'' 



dh'' 



logZ^'^^'i(Ext,/?0 



< [ci'^l+0(e)]|Ext|<[(Co + 0(e))|Ext|]l'=l (D.9) 



provided h e Uo. Combining (D.9) with the bound 

|Z^-^"(Ext,/lo)| < e-/r"'|Ext|gO(e)|aExt|g7riV,,(Ext)^ 

where /g™^'^ is the free energy of the contour model with activities 

K'^(Y) if y is small. 



(D.IO) 



if Y is large. 



(D.ll) 



we obtain 

d'' 



Z^"^i\Ext, h) 



dh'' 

On the other hand 

d'' 



< \k\l [{Co + 0(e))| Ext e-^r^"l^^^le^(^)l^^^^le^^^«-(^^^) 



h=hn 



(D.12) 



dh'' 



< |A;|!(Co|Xi|)'''le"(^~^('^)l^^le~-^l^*l<*+-^»^('"PP^*) < 
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Combining (D.7) with the inductive assumption (D.3) and the bounds (D.12) and 
(D.13), we may continue as in Appendix C to get 



<|ifc|!(Co+0(e))|l1'|)l'le--fl»'le^'''»''("')eOWI«"'lx 



h=ho 



X i: 



g-^|Ext| 



n 



(D.14) 



where now 



f = T — 6/e — 1. 



(D.15) 



Note the extra term 6/e with respect to (C.27), which comes from the term \k\/e in (C.13) 
(recaU that we assumed \k\ > 6). 

Given the bound (D.14), the proof of (D.3) for v{W) = n now foUows using Lemma C.2 
from Appendix C. This concludes the proof of Lemma D.2. □ 

Proof of Lemma 4.7 iv). 

Starting from (D.14), statement iv) of Lemma 4.7 is obtained in the same way as 
statement v) of Lemma 4.6 was obtained in Appendix C. □ 

As a coroUary of this proof, one obtains the analogue of (C.32) and (C.33) for derivatives, 
namely 



Z,{W,h) 



< \k\\ (Co + 0{e))\W\f\ e^^^«-(^)e(0(^)+0(^""'")l^^le-^l^l x 



h=ho 



X max e ^ 
ucw 



-^\W\U\^-{r/4Cs)\dU\ 



(D.16) 



and 



dh^ 



Z,{W,h) 



< \k\\ {Co + 0{e))\W\f\ e^rNev{W)^iOie)+Oie-^-/^))\dW\^-f\W\ ^ 



h—hr, 



xmax{e-^l^l,e-(2^-/4^3)|T^|''/(''-)| . (djj) 
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Appendix E: Proof of Lemma 5.1 and 5.2 

Proof of Lemma 5.1. 

Observing that all components W of Int'^^-' Ya obey the bound \W\ < max | the 

statement i) of Lemma 5.1 immediately follows from Lemma 4.6. 
In order to prove ii), we first note that for Ya — {Yi, • • • , Y^}, 



piyA)e 



Eg(supp Ya) _ 

n 



which implies that 

X g(e9-eo)(|suppyAU + |IiitYAnsuppA|d) ^ (E.2) 



Next we use Lemma 4.6 to bound 

N 



n 

m=l 



Z^ilntl^^YA.h) 



Z'(lntl^^YA,h) 



< gajlnt'o) yAH-0(e)|aint(°) yA|g27r7Va(Int(°) Ya) _ ^£^3^ 



Bounding Cq — eo < aq + 0(e), observing that 

I Int^°) Ya\ + \ Int Ya n supp A| = | Int Ya | , (E.4) 

and bounding |9Int'-°'' Ya\ < I^IntYAl + |9suppyl| < I^IntYAl + 2(i|suppA|, we get the 
bound 

Bounding now ATg (Int Ya) by Ci|Ya|, and observing that Hyer^ XqO^) implies that 
aq\V{YA)\ < (a + 1 + 0{€))\Ya\ due to the bound (4.26), we finally get 

|-K'g(>A)| < CAe"^^^^"^^+^^i^^)l^^le^^+"+^^')^l'^'*le^(')l"'PP^I , (E.6) 

which implies the bound (5.15). 

We are left with the proof of iii). By (E.l), Lemma D.l and the assumptions (3.9) and 
(3.25b), 



< \k\\CAC^o^ (IsuppYaI + |Ya| +2|IntY4n suppA])''"' 



X 



X g-(T-7T)|>A|g27TJVa(IntyAnsuppA)g(eq-eo)(| suppYaU + I IntYAHsuppAld) ^ 
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On the other hand, 



yr Z^(Intff yA,/l) 



< ((2Co + 0(e))| Int(°) Ya\) e^^^'^''"' x 



X e 



0(e)|aint(°) Ya L27TAra(Int(°) Ya) 



(E.8) 



by Lemma 4.7, while 



n ^^(^) 



YeVA 



/_ X |fe| 



(E.9) 



by (D.6). 

Combining the bounds (E.7) through (E.9), and bounding 
Co (Isupp YaI + + 2| Int Ya n supp A|) + (2Co + 0(e)) | Int(°> Ya\ + Cx J] [^(i^)! < 

<Co|supp^| + (2Co + 0(e))(|y^| + |Inty^|) + Ci in>^)|<|supp^|eOWI^^I, 

re^A (E.IO) 

we may then continue as in the proof of (E.6) to obtain the bound (5.16). □ 

Proof of Lemma 5.2. 

We start from the representation (5.4) and use the assumptions (3.25) and (3.11) to bound 



< C A^^'^^^^^'^'^^^^^ er'^'^~^'^^^'^^'^^^e^^'^'^\^'^^^ , (E.ll) 
Lemma 4.6 to bound 

N 



n Zm(IntS)r^,/i) 



m=l 



< -/|Int(")yA|gO(e)(|rA| + |suppA|) 7rJVa(Int^"JyA) 



(0) 



(E.12) 



and the inequahty (C.32) in conjunction with the estimate 

I^Ext^"^ YaI < l^suppAl + lay I + \dV{YA) \ < 2d\suppA\ + \dV\+C3\YA\ 
to bound 



|Z,(Ext(o)yA,/i)| <eO(^')(l^^l+l^"PP^I)e^"^«(^^*'°'^^)e-^l^^^'°'^^lx 



X e*^*^*^^!^-*' max e * 

U"CExt(o) Ya 



^1 Ext(o) YAW|g-(r/4C3)|a[/| 



, (E.13) 
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where e is the constant introduced in Lemma 5.2. Combining the bounds (E.ll) through 
(E.13) with (5.4), we obtain 



X 



C7CExt(0) Ya 



where we used the bounds aq\ Extn supp^| + (aq/4)| Ext^°) 1^ \ < iaq/4)\ ExtY^ \ U\ 
and Nq{V) < \dV\. Extracting the factor 

maxe-^l^-l max g"^! ^^*^-\^le-(^/^^3)|at/| 

Ya t/cExt(") Ya 

< maxe-(^/^^^3)|SV(i^A)| f,-^\{v\v{YAmu\^-{T/4Cs)\du\ 

Ya [/cExt(o) Ya 

< maxe-^l^\^le-(-/4^3)|9S| < ^^L-'^m ^-ir/4Cs)\dva 

~ SCV ~ I J 

from the right hand side of (E.14), we are left with a sum YIya e^*-^'^/^""*^-"^^' which we 
bound as follows 



^g-(3r/4-i)|r^| < ^ _L J2 e-(3T/4-l)|y| 

Y:V{Y)nsuppA^0 



n 

Ya n=0 



n 



< gO(?)| supp^l 



Putting everything together, we obtain the bound i) of Lemma 5.2. 

In order to prove ii), we generalize (E.ll) through (E.13) to derivatives. In (E.ll), these 
derivatives produce an extra factor, 

CA\k\\ (ColsupplOil + Col supp^ \ supp YaI)'''' < CaM (2Co|supp A| + Coir^l)'^' 

<C^|/c|!(Co|supp^|)l^le^WI^-l, 

while in (E.12) and (E.13), they produce factors 

Ifel 



|/c|!((Co + 0(6))|Int(°)y^|) 



and 



\k\\{{Co + 0{e))\Ext^^^YA\) . 
On the right hand side of (E.14), this leads to an extra factor 

C^|A;|!((Co + 0(e))|V|)l'=leOWI^-l. 

Observing that the sum g-^^^/^-^^^))!^^! can be bounded by e°(^)l^"PP^I as weU, we 
obtain Lemma 5.2 ii). □ 
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